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Development of Unsteady Aerodynamic Analyses for Turbomachinery 
Aeroelastic and Aeroacoustic Applications 


Summary 


Theoretical analyses and computer codes are being developed under Contract N ASS- 
25425 for predicting compressible unsteady inviscid and viscous flows through blade rows. 
Such analyses are needed to determine the impact of unsteady flow phenomena on the struc- 
tural durability and noise generation characteristics of turbomachinery blading. Emphasis 
is being placed on developing analyses based on asymptotic representations of unsteady flow 
phenomena. Thus, flows driven by small-amplitude unsteady excitations in which viscous 
effects are concentrated in thin layers are being considered. The resulting analyses should 
apply in many practical situations, lead to a better understanding of the relevant physics, 
and they will be efficient computationally, and therefore, appropriate for aeroelastic and 
aeroacoustic design applications. 

Under the present phase (Task III) of this contract research program, the effort has been 
focused on providing inviscid and viscid prediction capabilities for subsonic unsteady cas- 
cade flows. Thus, a composite (global/local) mesh, gust response solution capability for 
the linearized inviscid component of the flow has been developed, and an existing nonlin- 
ear viscous-layer analysis has been extended and coupled to the inviscid analysis to provide 
a weak inviscid/ viscid interaction solution capability for unsteady cascade flows. The lin- 
earized inviscid and the nonlinear viscous layer analyses are described in detail in this report. 
Numerical results are also presented to illustrate the effects of vortical excitation on the invis- 
cid pressure response and acoustic excitation on the viscous- layer response for representative 
two-dimensional compressor, turbine and flat plate cascade configurations. These results 
show that nonuniform mean flow phenomena can have a substantial impact on the response 
of cascades to vortical gusts, particularly for gusts occurring at high wave numbers. In 
addition, for flat plate cascades the viscous-layer responses to acoustic excitations from up- 
stream tend to be linear, but high frequency or high amplitude acoustic excitations from 
downstream, which travel against the main stream flow direction, are likely to excite signifi- 
cant nonlinear responses. Finally, because of boundary layer separation, weak inviscid/viscid 
interaction solutions can only be determined for a very restricted range of flow conditions. 
Therefore, a strong, or simultaneous, coupling of the inviscid and viscous-layer analyses must 
be developed to provide a useful unsteady viscous aerodynamic analysis for realistic cascade 
configurations. 
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1. Introduction 


The unsteady aerodynamic analyses intended for turbomachinery aeroelastic and aeroa- 
coustic predictions must be applicable over a wide range of operating conditions and to un- 
steady excitations at moderate through high reduced frequencies. In addition, these analyses 
must be capable of predicting unsteady pressure responses that arise from various sources of 
excitation. The latter include structural (blade) motions, variations in total temperature and 
total pressure (entropy and vorticity waves) at inlet and variations in static pressure (acous- 
tic waves) at inlet and exit. Finally, because of the large number of controlling parameters 
involved, there is a stringent requirement for computational efficiency, if the analysis is to be 
used successfully in the blade design process. To satisfy this latter requirement a number of 
restrictive assumptions must be introduced into the development of an appropriate unsteady 
aerodynamic model. 

The theoretical analyses that have been developed to predict the aeroelastic and aeroa- 
coustic behavior of turbomachinery blading, i.e., the onset of blade flutter, the amplitudes 
of forced blade vibration and the sound pressure levels upstream and downstream of the 
blade row, have, for the most part, been based on the following geometric and aerodynamic 
assumptions. The blades of an isolated, two-dimensional cascade are usually considered with 
the effects of neighboring structures being represented via prescribed nonuniform flow condi- 
tions at inlet and exit. In addition, the Reynolds number is assumed to be sufficiently high 
(i.e., Re — > oo) so that the displacement of the inviscid flow by viscous layers has a negligible 
impact on the unsteady pressure field. Finally, the unsteady excitations are assumed to be 
sufficiently small so that a linearized treatment of the unsteady inviscid flow is justified. 

Until recently, the inviscid unsteady aerodynamic analyses that have been available for 
turbomachinery aeroelastic applications were based on classical linearized theory, as reviewed 
by Whitehead in Ref. [l]. Because of the limitations in physical modeling associated with the 
classical linearization researchers have developed more general inviscid unsteady aerodynamic 
linearizations. These account for the effects of important design features such as real blade 
geometry, mean blade loading and operation at transonic Mach numbers on the unsteady 
aerodynamic response of two-dimensional cascades [2, 3, 4]. The unsteady disturbances 
are regarded as small-amplitude harmonic fluctuations relative to a nonuniform steady (in a 
coordinate frame attached to the blade row) background flow. The steady flow is determined 
as a solution of a nonlinear inviscid equation set, and the linearized unsteady flow, as the 
solution of a set of linear equations with variable coefficients that depend on the underlying 
steady flow. This type of analytical model has received considerable attention in recent years, 
and we refer the reader to the articles by Whitehead [5] and Verdon [6, 7] for a description 
of the theoretical formulation and recent numerical results. 

In particular, the linearized inviscid flow analysis (LINFLO), described in Refs. [6] and 
[7], employs a velocity decomposition introduced by Goldstein [8, 9], and later modified by 
Atassi and Grzedzinski [10], to split the linearized unsteady velocity into rotational and 
irrotational components. This decomposition leads to a very convenient description of the 
linearized unsteady perturbation — one in which closed form solutions can be determined for 
the entropy and rotational velocity fluctuations in terms of the drift and stream functions of 
the underlying steady flow and numerical field methods are required only to determine the 
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unsteady potential. The latter is governed by an inhomogeneous wave equation in which the 
source term depends only upon the rotational velocity field. Finite- difference procedures for 
solving the potential equation have been developed for the prediction of unsteady subsonic 
and transonic flows excited by blade vibrations or acoustic disturbances [11], and for the 
prediction of unsteady subsonic flows excited by entropic and vortical gusts [12]. 

The unsteady flows of practical interest usually occur at high, but finite Reynolds number, 
so that viscous effects can have an impact on the unsteady pressure response. Provided that 
large-scale flow separations do not occur, the overall flow field can be separated conceptually 
into “inner” viscous or dissipative regions, consisting of thin layers that lie along the blade 
surfaces (boundary layers) and extend downstream from the blade trailing edges (wakes), and 
an “outer” inviscid region. Solutions to the complete flow problem can then be determined 
by an iterative process involving successive solutions of the inviscid and viscid equations. 
If the inviscid/ viscid interaction is “weak”, then at each step of the iteration, the inviscid 
and viscid solutions can be determined sequentially with the pressure being determined by 
the inviscid flow. However, in most flows, strong inviscid/viscid interactions occur due, 
for example, to boundary-layer separations, shock /boundary-layer interactions and trailing- 
edge/near-wake interactions, and the pressure must be determined by solving the inviscid 
and viscous layer equations simultaneously. 

To date, there has been no effort to couple a viscous-layer analysis to the steady and 
linearized unsteady inviscid aerodynamic analyses that have been developed for turboma- 
chinery aeroelastic and aeroacoustic applications. As a step toward this goal, a fully nonlin- 
ear uns teady vis cous-layer analysis (UNSVIS), which can also be applied to predict steady 
flows, was developed in Ref. [13]. This approach, rather than one in which the results of 
separate nonlinear steady and linearized unsteady viscous analyses are superposed, will allow 
an assessment to be made of the relative importance of nonlinear unsteady effects in vis- 
cous regions. In the UNSVIS analysis the flow in the viscous layer is described by Prandtl’s 
equations, with algebraic models used to account for the effects of transition and turbulence. 
These equations are solved in terms of Levy- Lees type variables using a finite-difference tech- 
nique in which solutions are advanced in time and in the streamwise direction. Numerical 
solutions are determined by marching implicitly, first in time and then in the streamwise 
direction, over several periods of unsteady excitation from an initial steady solution, and 
from an approximate time- dependent, upstream flow solution. Numerical results were pre- 
sented in [13] for flat-plate airfoils and for the turbine cascades, studied experimentally by 
Dring et al. [14, 15]. These were found to be in good agreement with previous analytical and 
numerical solutions for flat-plate airfoils and with the experimental results for the turbine 
cascades. 

1.1 Scope of the Present Effort 

The objective of the research program being conducted under Contract NAS3-25425 is 
to provide efficient theoretical analyses for predicting compressible unsteady flows through 
two-dimensional blade rows. Such analyses are needed to understand the impact of unsteady 
aerodynamic phenomena on the aeroelastic and aeroacoustic performance of the blading. The 
work being conducted under this contract is directed primarily towards low-speed aeroelastic 
applications, however, for the most part it will apply more generally to the aeroelastic and 
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aeroacoustic performance of turbomachinery blading operating at high subsonic through 
transonic Mach numbers. 

In the first phase of this program [16] the linearized inviscid analysis (LINFLO) was 
extended to predict the response of a cascade to entropic and vortical excitations. Here, the 
velocity decomposition of Refs. [8] — [10] was employed to split the unsteady velocity into ro- 
tational and irrotational parts. Closed form solutions were then determined for the unsteady 
entropic and rotational velocity fluctuations, and the potential fluctuations were determined 
numerically on an H - type mesh in which the streamlines of the steady background flow were 
used as mesh lines. Numerical solutions were reported for several configurations including 
flat-plate cascades, a compressor exit guide vane, a high-speed compressor cascade, and a 
turbine cascade. 

Under the present effort the linearized inviscid analysis (LINFLO) and the nonlinear 
viscous-layer analysis (UNSVIS) have been extended and coupled to provide a weak vis- 
cid/inviscid interaction solution capability for unsteady cascade flows. In particular, a com- 
posite (global/local) mesh solution capability for flows driven by entropic and vortical gusts 
has been incorporated into LINFLO. The meshes used in the composite analysis are the gen- 
eral curvilinear H - mesh of Refs. [12] and [16] and a polar type local mesh. Also, analytical 
expressions for the values of the inviscid flow variables at moving blade and wake surfaces, 
i.e., at the edge of viscous layers, have been developed and incorporated into LINFLO, to 
provide the inviscid information needed for an unsteady viscous-layer calculation. 

The UNSVIS analysis has been extended so that the unsteady viscous effects in the 
vicinity of leading-edge stagnation points and in blade wakes can be predicted. In particular, 
a similarity analysis has been developed to predict unsteady viscous compressible flow in the 
vicinity of a moving leading-edge stagnation point and incorporated into the existing UNSVIS 
code. This analysis provides the “initial” upstream information needed to advance or march 
the viscous-layer calculation downstream along the blade and wake surfaces. Initial upstream 
profiles for the latter calculation are specified at a distance from the mean stagnation-point 
location which is beyond the furthest point at which any reverse flow, associated with the 
motion of the stagnation point, might occur. In addition, the wake analysis used previously 
in UNSVIS has been extended so that the changes or jumps in the inviscid velocity that 
occur across vortex-sheet unsteady wakes could be properly accommodated. 

The linearized inviscid and nonlinear viscous-layer unsteady aerodynamic analyses are 
described in detail in this report, and demonstrated via application to several cascade con- 
figurations. In particular, the inviscid analysis has been applied to predict the pressure 
response of compressor and turbine cascades, and their flat plate counterparts, to vortical 
excitations. Our numerical results indicate that it is important to account for mean flow 
nonuniformity in predicting the pressure response resulting from the interaction of a vorti- 
cal gust and a cascade, particularly for high-wave-number gusts. The nonlinear unsteady 
viscous-layer analysis has been applied to study the viscous-layer response of an unstaggered 
flat-plate cascade and a turbine cascade to external pressure excitations. The flat-plate re- 
sults indicate that the viscous-layer responses to upstream pressure excitations are essentially 
linear over a broad range of excitation amplitudes and temporal frequencies, but significant 
nonlinear responses occur for downstream pressure excitations at high amplitude or high 
temporal frequency. Finally, the coupled LINFLO/UNSVIS analysis has been applied to a 
turbine cascade subjected to a pressure excitation from upstream to demonstrate the current 
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weak inviscid/viscid interaction solution capability 


on a realistic cascade configuration. 
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2. Physical Problem 


We consider time-dependent, high Reynolds number (Re) flow, with negligible body 
forces, of a perfect gas with constant specific heats and constant Prandtl number (Pr) through 
a two-dimensional cascade, such as the one shown in Figure 1. The unsteady fluctuations in 
the flow arise from one or more of the following sources (Figure 2): blade motions, upstream 
entropic and vortical disturbances, and upstream and/or downstream acoustic disturbances 
which carry energy toward the blade row. These excitations are assumed to be of small 
amplitude, periodic in time, and periodic in the cascade “circumferential” direction. 

The mean or steady-state positions of the blade chord lines coincide with the line segments 
■q = £tanO + mG, 0 < £ < cos0 ,m = 0, ±1, ±2, ... , where £ and 77 are coordinates 
in the axial-flow and cascade directions, respectively, m is a blade number index, 0 is the 
cascade stagger angle, and G is the cascade gap vector which is directed along the 77- axis 
with magnitude equal to the blade spacing. The blade motions, '^-B m (x, t), where is 

the displacement of a point on a moving blade surface ( B m ) relative to its mean or steady- 
state position (B m ), x is a position vector and t is time, are prescribed functions of x and 
t. The entropic, S_oo(x,<), vortical, ^.^(x, f), and acoustic, Pi,?oo(x,t), excitations, where 
the subscripts —00 and +00 refer to the far upstream and far downstream flow regions, 
respectively, are also prescribed functions of x and t. However, a prescribed far-field behavior 
must be a solution of the fluid-dynamic field equations, and the incident pressure disturbances 
must satisfy the requirement that acoustic energy travels toward the blade row. Since the 
unsteady aerodynamic excitations are periodic in 77 and t , it is often useful to describe 
them via Fourier series representations. In this case the complex amplitudes, fundamental 
frequencies and fundamental circumferential wave numbers of the various excitations are the 
prescribed quantities. 

We assume that far upstream (say £ < £_) and far downstream (£ > £+) from the 
blade row, the mean flow in the inviscid region is at most a small irrotational perturbation 
from a uniform free stream. In addition, any shocks that might occur are assumed to be 
of weak to moderate strength. Finally, blade shape and orientation relative to the inlet 
freestream direction, the inlet to exit mean static pressure ratio and the amplitudes, modes, 
frequencies and wave numbers of the unsteady excitations are such that viscous effects are 
confined within thin layers, which remain essentially attached to the blade surfaces, and 
extend downstream from the blade trailing edges. 

In the present discussion all physical variables are dimensionless. Lengths are scaled 
with respect to blade chord, time with respect to the ratio of blade chord to upstream 
freestream flow speed, density and viscosity with respect to their upstream freestream values, 
respectively, velocity with respect to the upstream freestream flow speed, normal and shear 
stresses with respect to the product of the upstream freestream density and the square of 
the upstream freestream speed, temperature with respect to the square of the upstream 
freestream speed divided by the fluid specific heat at constant pressure, and entropy with 
respect to the fluid specific heat at constant pressure. 
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2.1 Mathematical Model 


The field equations that govern the fluid motion follow from the mass, momentum and 
energy conservation laws and the thermodynamic relations for a perfect gas. After ensemble 
averaging the resulting equations, equations for the (statistical) mean values of the flow 
variables are determined. These contain turbulent correlations, which must be estimated 
empirically. As a consequence of the high Reynolds number assumption, separate sets of 
equations can be constructed to describe the flow in the “outer” inviscid region and in the 
“inner” viscous- layer regions. The flow in the inviscid region is governed by a coupled set 
of nonlinear differential equations (i.e., the Euler equations) along with jump conditions at 
shocks and at vortex-sheet unsteady wakes. In continuous subregions the energy equation 
can be replaced by the requirement that the entropy following a fluid particle must remain 
constant. In addition, the attached flow assumption requires that the inviscid unsteady flow 
must be tangential to the moving blade surfaces, and the uniform flow conditions at inlet 
and exit, the entropic and vortical fluctuations at inlet, and the static pressure disturbances 
at inlet and exit that carry energy towards the blade row must be specified. The flows in 
the viscous layers are governed by Prandtl’s equations and are subject to no-slip conditions 
at the moving blade surfaces. In addition, the streamwise velocity and the thermodynamic 
properties of the fluid at the edges of the viscous layers must be matched to the corresponding 
inviscid quantities at the blade surfaces and along the reference wake streamline. 

Even with the simplifications resulting from the high Reynolds number approximation, 
the unsteady aerodynamic problem still involves prohibitive computational times. Because 
of this, the traditional approach has been to examine linearized forms of the governing invis- 
cid equations with the intention of providing useful theoretical predictions of the unsteady 
aerodynamic response information required for turbomachinery aeroelastic and aeroacoustic 
applications. In one such approach, the unsteady flow in the inviscid region is regarded 
as a small perturbation of a nonlinear mean or steady background flow. In this case the 
first-order unsteady fluid motions are governed by linear equations with variable coefficients 
that depend upon the steady background flow. In addition, the first-order motions that arise 
from the various Fourier modes of unsteady excitation are not coupled. Hence, they can be 
determined independently. Indeed, it is sufficient to develop solution procedures for a single 
harmonic (in t and 77) component of a given disturbance. Solutions for arbitrary disturbances 
and arbitrary combinations of various disturbances can then be obtained by superposition. 

Thus, in linearized inviscid analyses we consider small-amplitude blade motions of the 
form 

”^B m (x + mG, 1) = Re{r B (x) exp[i(wi + m<r))}, x e B, (2.1) 

where rjj is the complex amplitude of the blade displacement, a is the phase angle between 
the motions of adjacent blades, Re{ } denotes the real part of { } and B denotes the reference 
(m = 0) blade surface. The external aerodynamic excitations are small-amplitude, harmonic 
perturbations of a uniform freestream flow, and therefore must be of the form 

s(x, t ) = S_oo(x - V.oot) = Refs-oo exp[i(K_oo • x + Ujt)]} . . . , ( < , (2.2) 

COM) = C-oo( x “ V -M) = IMC -00 ex P[*(*-°° -x + wf)]}... , (2.3) 
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and 


p/(x, t) = pi iT oo(x, t ) = Re{p/, Too exp[-^ Too <f + *’(#c Too • x + u;*)]} . . . , £ > 6 f • (2.4) 

Here s_ <», C-«> and p/ iToo are the complex amplitudes of the entropic, vortical and acoustic 
excitations, u is the temporal frequency, #c Too is the wave number with oo = \ 

and the /3 Too are attenuation constants. It follows from the governing equations that the 
temporal frequency and wave number of an entropic or vortical excitation are related by 
u — _ k_oo • V-oo, where 'V-oo is the uniform relative inlet velocity, but more complicated 
relationships exist between u) and k Too for pressure excitations [17]. 

As a consequence of the foregoing unsteady aerodynamic linearization, the first-order 
inviscid flows arising from excitations of the form (2.1) through (2.4) will be harmonic in 
time, and they will satisfy a blade-to-blade periodicity condition. Therefore, these flows 
can be determined by solving a time-independent, linear, boundary-value problem, over a 
single extended blade-passage solution domain, for the complex amplitudes of the linearized 
unsteady flow variables. Also, because of our assumptions regarding shocks and the steady 
flow far upstream of the blade row, the steady background flow will be isentropic and irrota- 
tional. In this case, closed form solutions can be determined for the unsteady entropic and 
vortical perturbations. Therefore, numerical field methods are required only for determining 
the first-order unsteady pressure fluctuations. 
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3. The Flow in the Inviscid Region 


3.1 Unsteady Perturbations of a Potential Mean Flow 


For small amplitude unsteady excitations, i.e., |s_oo|, IC-ooli e ^ c -5 ~ C*( e ) << 1* the 
time- dependent flow can be regarded as a small perturbation of an underlying nonlinear 
mean or steady background flow. Thus, for example, we can express the time-dependent 
velocity in the form 

V(x,f) = V(x) + v(x, t) + ... , (3.1) 

where V(x) is the local mean velocity and v(x, t) is the first-order (in e) unsteady velocity. 
In addition, Taylor series expansions, e.g., 

V S = (V + K-VV + ...)|s , (3.2) 


and surface-vector relations, i.e., 


ts = t s + 



n s + ... 


and ns = ns 



(3.3) 


can be applied to refer information at a moving blade, wake or shock surface to the corre- 
sponding mean surface. Here 7?.(x, t) = xs — x s is the displacement of a point on the moving 
surface relative to its mean position, the subscripts S and S refer to the moving and mean 
surface locations, respectively, and the unit tangent, T, and normal, n, vectors are oriented 
such that n x r points out from the page 

The equations that govern the steady and first-order unsteady flow properties are de- 
termined by substituting the foregoing expansions into the full time-dependent, nonlinear 
governing equations, collecting terms of like power in e and neglecting terms of 0(e 2 ) or 
higher. This procedure provides a nonlinear boundary-value problem for a steady back- 
ground flow and a linear variable-coefficient problem for the first-order unsteady flow in 
which the coefficients depend on the underlying steady flow. The first-order or linearized 
unsteady fluctuations caused by a periodic harmonic unsteady excitation will be harmonic 
in time, e.g., 

v(x, t) = Re{v(x)exp(tu>f)} . (3.4) 

In addition, the steady and, for an excitation at wave number k ,, = crG -1 , the linearized 
unsteady flow properties will satisfy the following circumferential conditions, e.g., 


V( x + mG) = V(x) and v(x + mG) = v(x) exp(imo-) . (3.5) 


Thus, solutions to time-independent nonlinear steady and linearized unsteady flow problems 
are required only over a single extended blade-passage region of the cascade. In addition, 
since analytic far-field solutions can be determined [17], the numerical solution domain can 
be restricted further to a single extended blade-passage region of finite extent in the axial 
direction. 
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The Steady Background Flow 

Because of our assumptions regarding shocks and the flow far upstream of the blade row, 
the steady background flow is isentropic and irrotational; i.e., V = V$, where $ is the 
steady velocity potential. Thus, the field equations that govern this flow are 

V • (pV$) = 0 (3.6) 


and 


(M-ocV/M ) 2 = (M.ooA) 2 = p^ = (7 = (7 - m 2 -ooT 


= 1 - 


1 Ijv/ 2 f(V$) 2 - 1] = - 

2 M -ooU vv l M 2 + (7 — 1)M 2 


(3.7) 


where M , A, p, P and T are the local Mach number, speed of sound propagation, density, 
pressure and temperature, respectively, in the mean or steady background flow and 7 is the 
specific heat ratio of the fluid. Surface conditions for the zeroth-order or steady flow apply 
at the mean positions, B m , W m and Sh mi „, of the blade, wake and shock surfaces, where the 
subscript n refers to the nth shock associated with the ntth blade. Since, by assumption, 
the flow remains attached to the blade surfaces, a flow tangency condition applies at such 
surfaces. In addition, mass and tangential momentum must be conserved across shocks, and 


the steady pressure and normal velocity component must be continuous across blade wakes. 

Numerical procedures for determining two-dimensional steady potential flows through 
cascades have been developed extensively, e.g., see [3, 18], particularly for flows with subsonic 
relative inlet and exit Mach numbers (i.e., Af T<X) < 1). In such calculations far-field boundary 
conditions are imposed at axial stations placed at finite distances upstream and downstream 
(i.e., at £ = £j:) from the blade row, where linearized solutions describing the behavior 
of the steady potential can be matched to a nonlinear near- field solution. In addition, a 
Kutta condition is usually imposed at blade trailing edges in lieu of prescribing an exit 
freestream property. Finally, the usual practice is to solve the conservative form of the mass- 
balance equation (3.6) throughout the entire fluid domain. Thus, the shock- and wake-jump 
conditions are not imposed explicitly. Instead, shock phenomena are captured through the 
use of special differencing techniques. W&ke conditions are satisfied implicitly, because the 
fluid properties are continuous and differentiable across wakes. If needed, mean shock and 
wake locations can be determined a posteriori from the resulting steady solution. 


3.2 Linearized Unsteady Equations 

The field equations that govern the first-order unsteady perturbation of an isentropic and 
irrotational steady flow can be expressed as a system of coupled differential equations for the 
complex amplitudes of the first-order entropy (5), velocity (v) and pressure (p), respectively 
(see [8, 19]). In general, we require a solution to this system subject to prescribed flow 
conditions far upstream and far downstream from the blade row, the flow tangency condition 
at moving blade surfaces, and, since the inviscid field equations apply only in continuous 
regions of the flow, jump conditions at moving shocks and blade wakes. 

As indicated by Goldstein [8, 9] the system of field equations that govern the linearized 
unsteady flow can be recast into a very convenient form by decomposing the unsteady velocity 
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into rotational (v fi ) and irrotational (V^) parts. Under this decomposition the rotational 
velocity, v/<, is divergence- free far upstream of the blade row, i.e., V • Vr = 0 for £ < £-, and 
the unsteady pressure depends only upon the potential <j> through the relation p = — pD( j>/ Dt , 
where DjDt = iu + ■ V is a convective derivative operator based on the mean flow 

velocity. The unsteady vorticity is given by ( = V X v fi . Here we apply a modified form 
of the Goldstein decomposition, which was introduced by Atassi and Grzedzinski [10], to 
facilitate the numerical resolution of the velocity potential, i.e., we set 

v = v* + V<f> = Vfl + , (3.8) 

where <f>+ is a convected or pressure-less potential ( D<j> m /Dt = 0), which renders v, • n = 
(vr + V<^„) • n = 0 at blade and wake mean positions. 

With the decomposition (3.8) the system of field equations that governs the unsteady 
flow variables s, vr and <j) can be written in the form 

= 0 (3.9) 

Dt 

§-(v R - sV$/2) + [(v R - sV$/ 2) • V]V* = 0 (3.10) 

and 

-^(A- 2 ^) - p - 1 V • (pV<f>) = p • [p(y R + V^.)] = r'V ■ (Pv,) . (3.11) 

These equations are coupled only sequentially, i.e., they can be solved in order to determine 
the complex amplitudes s , v R and <f>. Moreover, since the entropy fluctuation is governed 
by a convection equation, it can be determined as a solution of Eq. (3.9) in terms of the 
prescribed upstream entropy field. The rotational velocity fluctuation, is governed by 
a modified convection equation and can be determined as a solution of Eq. (3.10) in terms 
of the prescribed upstream entropy and rotational velocity fields. The velocity potential, 
<j>, is governed by a second-order partial differential equation with source term p 1 V • (pv„) 
and depends, therefore, upon the rotational velocity, the convected potential, and boundary 
condition information both fax upstream and far downstream of the blade row as well as on 
the blade, wake and shock surfaces. 


Far-Field Behavior 

We have assumed that the mean or steady flow is at most a small (i.e., of O(e)) pertur- 
bation from a uniform stream both far upstream (£ < £_) and far downstream (£ > £+) from 
the blade row. Therefore, in these regions, the first-order (in e) unsteady field equations can 
be reduced to constant coefficient equations for which analytical solutions can be determined 
(see [17]). For example, it follows after replacing by V_oo in (3.9) and (3.10) that the 
complex amplitudes of the entropy and rotational velocity fluctuations far upstream of the 
blade row are given by 

•s(x) = s.oo exp (i#c_oo • x), £ < i- , (3-12) 
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and 


Vfi(x) = Vr-oo exp(i/c_oo • x), £ < , (3.13) 

where S-oo and = (V x v fl )_ ^ = ik_oo x \r-oo, the complex amplitudes (at x = 0) of 
the incident entropic and vortical gusts, are prescribed. The circumferential wave number, 
k»7,too, is <r/G and it follows from (3.9) and (3.10) that k • V_oo = — w. Since the rotational 
velocity is divergence-free far upstream of the blade row, the vectors #c_oo and Vj? t _ oo are 
orthogonal, and therefore, 

V« oo = *(*-«> x C-cc)/l*-oo| 2 • (3.14) 

The velocity potential fluctuations in the far upstream and far downstream regions de- 
pend upon the acoustic excitation as well as the response of the cascade to the imposed 
unsteady excitations. We can set 

4>{x) = 4 >e{x) + <M X ) for £>£,:, (3.15) 

where the potential component 4e accounts for acoustic excitations, i.e., pressure distur- 
bances that either attenuate as they approach the blade row or propagate and carry energy 
towards the blade row. In particular, it follows from (3.11) with = V Too that for an 
acoustic excitation at temporal frequency u>, 

<M X ) = o exp[/? Too ^ + *#c T oo • x], i > , (3.16) 

where 

4*1,^ oo = P^o {$foo C° s 4" ' V^oo)]} Pl&oo > (3-17) 

the complex amplitudes Pi^oo are prescribed, «^ iToo = <r/G, and /? Too and depend 

upon the inlet/exit freestream conditions, and the temporal frequency (ur) and interblade 
phase angle ( a ) of the acoustic excitation. The potential component <f>R is associated with 
the acoustic and vortical response of the blade row and therefore, must be determined as 
part of the unsteady solution. 


3.3 Entropy and Rotational Velocity 


Closed form solutions for the linearized entropy and rotational velocity fluctuations can 
be determined in terms of independent variables that describe the steady background flow [8, 
12]. For this purpose we introduce the Lagrangian coordinate vector 

X = Ae r + tfe* . (3.18) 


Here, 


and 


A(x) = x_ + r V-'dT • , 

V x _+[«(x)-*(x_)]e w 

'P(x) = x_ + [ p(e z x V) • dr 
Jx-+e N 


(3.19) 

(3.20) 


are the drift and stream functions, respectively, of the steady background flow. Also, the 
unit vectors ey, ejy(= e 2 x ej) and e 2 point in the inlet freestream direction, normal to 
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the inlet freestream direction, and out from the page, respectively, x_ is the position vector 
to the point of intersection (£_,/?_) of the reference blade stagnation streamline and the 
axial line £ = dr* is a differential element of arc length along a streamline, and dr is a 
differential vector tangent to the path of integration. 

Note that DX/Dt = (V-V)X = = e T and X -► x, as £ -»■ -oo. Thus, the solution 

to the entropy transport equation (3.9), which satisfies the upstream condition (3.12), is 

s(x) = s_oo exp(f#c_ 00 • X) , (3.21) 

and the solution to the rotational velocity transport equation (3.10) , which satisfies the 
condition (3.13), is 

v H (x) = [V(X • A_ oo) + s-oo V$/2] exp(i»c_ <x> ■ X) , (3.22) 


where 

•A ~ co = V H,— oo s — <x>V— oo/2 • 

The complex amplitude of the unsteady vorticity is 

C = V x v R = V(iK _ 0 o • X) x [V( A_oo • X) + s_ooV$/2] exp(i/c_oo • X) , 


(3.23) 

(3.24) 


and C — *■ C-oo = iK x Vr,_ oo as X x. 

If the steady background flow stagnates at blade leading edges, as it will for realistic 
configurations, the drift function will have a logarithmic singularity at the mean blade and 
wake surfaces, i.e., A -* a 0 In n as n — ► 0, where n is the normal distance from the surface 
and a 0 is a real constant. As a result, the exponential function exp(fic_ 00 • X) will be 
indeterminate, and the normal component of the rotational velocity will be singular, i.e., 
Vh . n _+ ai n~ l exp(iK_oo • X), where a\ is a complex constant, at such surfaces. We can 
remove this singular behavior from the blade and wake surface conditions that determine 
the unsteady potential, <f>, by prescribing a convected potential of the form [10] 


<f>, = [-iu; -1 A.,*, • V-oo + F(V)] exp(i/c_oo • X) , (3.25) 


where 


*w = 


(/C-oo x A-oo) • e z Gcosf)- c 
2x(l — ia 0 ui) 


sin 


2?r['l'(x) — 'J'(x-)] 
G COS fl_oo 


(3.26) 


is a complex function that depends upon, among other things, the behavior of the mean flow 
in the vicinity of a leading-edge stagnation point. This choice of <j > , ensures that v* • n = 
(vr + V ^) • n = 0 at blade and wake mean positions. 

After combining (3.22), (3.25) and (3.26), we find that the complex amplitude of the 
source-term velocity, v, = v# + V </>„, is given by 


v. = 


'dF 


FV(tK_ oo . X) + ( — - w -1 (#6_ oo X A- oo) • e z ) + 5-00 V$/2 


) 


exp(f/c_ 0O • X) . 

(3.27) 

Note that v* behaves like s_ooV$ exp(t»c_oo • X)/2 in the immediate vicinity of the mean 
blade and wake surfaces, i.e., as n — * 0. Thus, v» • n = 0, but, if s_oo 0, the tangential 
component of the source-term velocity will be indeterminate at such surfaces. 
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The velocities v# and v* depend upon A and ^ and the first partial derivative of these 
functions. Therefore, the complex amplitudes of the unsteady vorticity, ( = V x = 
V x v*, and the source term, p _1 V • (pv*), in (3.11) depend also upon the second partial 
derivatives of A and Thus, an accurate solution for the nonlinear steady background 
flow is critical for determining the unsteady effects associated with entropic and vortical 
excitations. 

Numerical Evaluation 

The complex amplitudes of the entropy, rotational velocity, vorticity, and source term 
velocity are readily determined once the values of the drift and stream functions and their 
spatial derivatives are specified over the single extended blade-passage solution domain. 
For this purpose it is convenient to use an H - grid in which one set of mesh lines are the 
streamlines of the steady background flow. An H- grid which covers the solution domain, i.e., 
which is bounded by the upstream and downstream axial lines £ = ^ and two neighboring 
mean-flow stagnation streamlines, is appropriate. The locations of the latter are determined 
a posteriori from the solution for the nonlinear steady background flow. Once the boundaries 
of the i/-grid established, the locations of the interior grid points can be determined using 
an elliptic grid generation technique as described in Ref. [12]. 

Because a streamline mesh is used, the drift function can be evaluated at each point in 
the computational domain by a straightforward numerical integration of Eq. (3.19). The 
procedure used in Ref. [12] is simply to specify the drift function along the far upstream 
boundary £ = and then to evaluate this function along each streamline using a second- 
order accurate difference approximation. The derivatives of the drift and stream functions at 
a given grid point are determined using the finite difference operators developed by Caspar 
and Verdon [20]. Because the drift function is singular at blade and wake surfaces, one-sided 
difference approximations are used to evaluate the derivatives of this function at points on 
the mesh streamlines adjacent to these surfaces. 

3.4 Velocity Potential 

The unsteady potential ((f)) is determined as a solution of the field equation (3.11) subject 
to conditions at the mean blade and wake surfaces and in the far field. Flow tangency applies 
at the blade surfaces, the fluid pressure and normal velocity must be continuous across 
blade wakes, W m , and mass and tangential momentum must be conserved across shocks, 
Sh m}Tl . As a consequence of the small unsteady-disturbance approximation, the conditions 
on the linearized unsteady perturbation at the moving blade, shock and wake surfaces can be 
imposed at the mean positions of these surfaces (see Ref. [19]), with the mean wake (kF m ), 
i.e., the downstream stagnation streamlines, and shock (Sh mjTl ) locations being determined 
from the nonlinear steady solution. Thus, the flow tangency condition can be written as 

v • n = 'Vcf) * n = [ico r + (V<1> • t)(t * V)r — (r * V)V$] • n, x G B m (3.28) 

In addition, since the irrotational steady velocity and pressure are continuous and have 
continuous derivatives across the mean-flow downstream stagnation streamlines, the wake 
conditions reduce to 

[p] = {D(j)/Dt] = 0 and [v • n] = [V</>] • n = 0, x e W m , (3.29) 
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where the symbol [ J indicates the jump or change in a quantity at a surface at which the 
flow variables are discontinuous. Finally, if we neglect changes in entropy and rotational 
velocity across shocks, the conservation laws for mass and tangential momentum yield the 
following linearized shock-jump condition for a shock that terminates in the fluid 


I/W + Vfi - A-^ V$)] • n = -[/»][»w + (V* • t)t • V][([V*] • n)" 1 ^] ] 
- (V* • n)" 1 ^] r • V([p]V$ • r) , x 6 Sh m , n . 


(3.30) 


Equation (3.30) provides a relation for determining the jump in the unsteady potential, [<^], 
at the mean position of a shock. The shock displacement normal to the mean shock locus is 
then given by rsh • n = -([V$] • n)* 1 

The velocity potential in the far field is given by (3.15); the potential due to an acoustic 
excitation at frequency u> and circumferential wave number = 0*/G, by (3.16). Usually 

only acoustic excitations that are of propagating type are considered. For subsonic inlet and 
exit conditions ( M Too < 1) the velocity potential corresponding to a propagating acoustic 
excitation at temporal frequency u > 0 has the form (3.16) with /? TOO = 0, and 


(-) 


^TJ^OO ^ ^T)^00 


-1 


< K (+) 


(3.31) 


where the circumferential wave numbers at which cut-off occurs, /cj^oo, are given by 


= ‘*V-iM FOO (l - oo sin ± A - cos’ U T „) . (3.32) 


The axial wave number of the propagating acoustic excitation is 

oo) > (3.33) 

A 

where 

l^ool = 1(1 - cos’Ofoo)- 1 *; - 2 , (3.34) 

and 

= (wV^ + sinfl Too )/(l - M 2 ^ cos 2 H Too ) . (3.35) 

Analytic solutions to (3.11), with = V TOO , for the far-field potential component <f>fi 
which satisfy the requirements that acoustic response disturbances either attenuate with 
increasing axial distance from the blade row or propagate carrying energy away from or 
parallel to the blade row, and that vorticity must be conveeted downstream are given in 
Ref. [17]. These solutions contain arbitrary constants that are determined by matching the 
far-field analytic solutions for the velocity potential to near field numerical solutions. 


Numerical Evaluation 


A numerical resolution of the foregoing linear, variable-coefficient, boundary-value prob- 
lem for <f> is required over a single, extended, blade-passage region of finite extent. The field 
equation, Eq. (3.11), must be solved in continuous regions of the flow subject to surface 
conditions, Eqs. (3.28-3.30), at the mean blade, wake and shock surfaces. In addition, the 
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near-field numerical solution for the potential must be matched to far-field analytical solu- 
tions at finite distances (£ = upstream and downstream from the blade row. Numerical 
methods for determining <f> for subsonic and transonic flows have been reported in Refs. [2], 
[11], [20] and [21]. To date, transonic solutions have been determined only for flows in which 
only a single normal shock occurs in each blade passage and for unsteady flows in which 
s = £ = 0. Thus, numerical solution procedures, see [12] and [16], for the entropic and 
vortical gust problems have only been developed and implemented for subsonic flows. The 
development of such procedures for the transonic gust response problem remains, therefore, 
as a subject for future research. 

Because of the stringent and conflicting requirements placed on computational meshes 
for cascade flows, a composite- mesh, which is constructed by overlaying a polar-type local 
mesh and an H - type cascade mesh, has been adopted for determining the unsteady potential. 
The H mesh is used to resolve unsteady phenomena over the entire solution domain; the 
local surface-fitted mesh, to resolve phenomena in the vicinities of a rounded blade leading 
edge and/or a shock. The cascade mesh facilitates the imposition of the cascade periodicity 
conditions (3.5) and the matching of the analytic and numerical unsteady solutions at the far 
upstream (£ = £_) and far downstream (£ = £+) boundaries. Use of this mesh alone is often 
sufficient for resolving unsteady subsonic flows. The local mesh allows an accurate modeling 
of unsteady leading-edge and shock phenomena. It is constructed so that two “radial” lines 
coincide with the predicted mean shock locus to provide upstream and downstream shock 
mesh lines for the accurate imposition of shock-jump conditions. 

Since the cascade and local body-fitted meshes differ topologically, a zonal solution pro- 
cedure for overlapping meshes has been adopted in Ref. [11] for determining the unsteady 
potential. In the region of intersection between the two meshes, i.e., the region covered by 
the local mesh, certain cascade mesh points are eliminated depending upon their location 
within the local mesh domain. The discrete equations are written separately for the cascade 
and local meshes and coupled implicitly through special interface conditions, resulting in a 
single composite system of finite-difference equations that describe the unsteady flow over 
the entire solution domain. 

The finite-difference model used to approximate the unsteady equations on the cascade 
and local meshes has been described in detail in Ref. [20]. Algebraic approximations to the 
various linear operators, which make up the unsteady boundary-value problem, are obtained 
using an implicit, least- squares, interpolation procedure that is applicable on arbitrary grids. 
This procedure employs a nine point "centered” difference star at subsonic field points, 
and a twelve point difference star at supersonic points. At a blade boundary point a nine 
point one-sided difference star is used on the cascade mesh, while nine- or six-point one- 
sided stars are used on the local mesh. Shocks are fitted in the local-mesh calculation 
by approximating the shock-jump condition (3.30) using one-sided difference expressions to 
evaluate the normal derivatives of the unsteady potential on the upstream (supersonic) and 
downstream (subsonic) sides of the shock. At those points along the shock mesh lines at 
which the steady flow is continuous (i.e., at points lying beyond the end of the shock), the 
condition [<£] = 0 is imposed. 

The systems of linear algebraic equations that approximate the unsteady boundary- value 
problem on the cascade and local meshes are block-tridiagonal for subsonic flow and block- 
pentadiagonal for transonic flow. A subsonic solution on the H - mesh alone is determined 
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using a direct block inversion scheme. Composite (cascade/local) mesh subsonic or tran- 
sonic solutions are determined using a different scheme. Because of the cascade/ local mesh 
coupling conditions, the composite system of discrete equations contains a sparse coefficient 
matrix of large band width. Consequently, special storage and inversion techniques must be 
applied to achieve an efficient solution. Once the composite system of unsteady equations is 
cast into an appropriate format, it can be solved using Gaussian elimination [ 11 ]. 

3.5 The Inviscid Response 

At this point we have presented a linearized unsteady aerodynamic formulation that de- 
scribes the general first-order fluid-dynamic perturbation of an isentropic and irrotational 
subsonic or discontinuous transonic steady background flow. We have also outlined the nu- 
merical procedures used to evaluate the unsteady entropy, rotational velocity and velocity 
potential. Solutions to the linearized unsteady problem are required to determine the aero- 
dynamic response information needed for aeroacoustic and aeroelastic predictions, i.e., the 
unsteady pressure field and the unsteady pressures acting at the moving blade surfaces. In 
particular, we refer the reader to Refs. [7, 22] for a description of the various local and global 
unsteady aerodynamic response parameters that are used in aeroelastic investigations. 

Approximate solutions for the full time-dependent flow properties are constructed by 
superposing the results for the steady and the unsteady flow properties, e.g., 

^(X,0 = P(x)+p(x,t) + ')TPshm,n( X ’ t ) + -’- 1 ( 336 ) 

where p(x) = Re{p(x) exp(iu;t)}. The first two terms on the right-hand-side of Eq. (3.36) 
account for the steady and the first-harmonic contributions to the time-dependent fluid pres- 
sure acting at the point x. The third term accounts for the anharmonic contribution to the 
time-dependent pressure caused by the motions of shocks, and is determined by analytically 
continuing the steady solution from the mean to the instantaneous shock locations [22, 23]. 
Thus, the first-order anharmonic response depends on the steady values of the fluid proper- 
ties and the unsteady shock displacements normal to the mean shock loci, (72. • n)sh m ,„- The 
regions of anharmonicity are confined to thin strips containing the mean shock loci. 

The steady background flow is described by the velocity potential, $, which is determined 
in terms of a prescribed inlet Mach number (M-oo), flow angle (Q_oo), and cascade geometry. 
The steady velocity, V = V$, is determined from this potential, as are the steady values of 
the thermodynamic properties of the fluid [cf. (3.7)]. The total enthalpy, H = T + V 2 /2, of 
the steady background flow is constant and is given by 

H = = (7 - l)-'ACi + \ (3.37) 

The linearized unsteady flow is described by the dependent variables s, v„, <f> and • n, 
which are determined in terms of prescribed values of the frequency (u>), interblade phase 
angle (er), and the complex amplitudes, rg, S-,*,, and p/,joo» of the imposed unsteady 
excitation. The first-harmonic unsteady velocity is v = v„ + V<^, and the first-harmonic ther- 
modynamic properties of the fluid can be determined in terms of the independent variables 
5 and <f>. 
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It follows from the equation of state for a perfect gas and the fundamental thermodynamic 
identity relating entropy, pressure and density that 

P = 7 -1 (7 - l)(pt + pT) i ( 3 - 38 ) 

and 

s = r'p/p-p/p, ( 3 - 39 ) 

respectively, where t is the linearized unsteady component of the fluid temperature. After 
combining (3.38) and (3.39) and recalling that p = —pD<j>/Dt = —p(d/dt + • V)^, we 

find that the first-harmonic flow properties are related by 

pip + 5 = 7 ~ l p/P = ( 7 - 1 )-\ifT -s) = —A~ 2 D4>/Dt . (3.40) 

The first-harmonic total enthalpv h , is given by 

h — Ts Dt V $ • [v* -(- V <^>] (3.41) 

The Inviscid Flow Along Moving Blade and Wake Surfaces 


The values of the inviscid flow quantities along the moving blade and wake surfaces must 
be available to perform an analysis of the unsteady flows in the viscous layers. These can be 
determined in terms of their values at the mean surface locations and the blade and wake 
displacements via Taylor series expansions, e.g., 

Ps = Ps+Ps + K-VP\ s + 'i ! rPsh n + ... , ( 3 - 42 ) 


where 72 is the surface displacement vector. In addition, relations between the unit tangent 
and normal vectors at the instantaneous and mean surface positions, cf. (3.3), are needed to 
determine the velocity components at the moving surface from information available at the 
mean surface. Here, the subscript <5 refers to points Xs = x$ + 72-s on the moving surface, 
and S, to the points Xs on the mean surface. 

If we now restrict our consideration to continuous, i.e., shock-free steady background 
flows, the unsteady velocity at a moving blade or wake surface is given by 


V 5 = [V + v + (72 • V)V] S + ... 

= [V$ + v, + + (72 • V)V$] S + ... 


(3.43) 


Expressions for the tangential and normal components of the fluid velocity at a moving blade 
and wake surface follow from (3.3), (3.43) and the conditions • n$ = 0 and v, • ng = 0. 
We find that 

\ s T S = [V$ + v, + V<£ -f (72 • V)V$] S • r s -f ... ( 3 -44) 

and qk 

\s • ns = [V</> + (72 • V) V $]s • ns + ... = • ns . (3.45) 

The thermodynamic properties at the moving blade and wake surfaces are determined 
using Taylor series expansions along with .the thermodynamic property relations for the 
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steady (3.7) and linearized unsteady (3.40) flows. After performing the necessary algebra, 
we find that for shock-free flows 

i + 5 - s = ^ = (7-1)-' - «) = -V (§7 +'*■ V[(V*)V2]) s + ... . (3.46) 

Finally, the total enthalpy of the fluid at a moving blade or wake surface is given by 

H s = [H + h + K- VH] S + ... • (3-47) 


Since the total enthalpy is constant in the steady background flow, i.e., H — -ff-oo, an< f 
[V$ • v.] s = 5(V$ • r) 2 / 2 at blade and wake mean positions, we can write 


H s = H- 00 + T s s s -[D4>/Dt]s + V<f>-(V4> + v.)s+ ••• 
= [77-00(1 + s) - Dtp/Dt + • V^] s + ■ • • • 


(3.48) 


Expressions for the anharmonic values of the time-dependent flow properties in the field 
and at moving blade and wake surfaces can be readily determined following the methods 
of Refs. [22] and [23]. However, in the present effort only subsonic inviscid and viscous- 
layer flows will be calculated, and therefore, only the steady and first-harmonic values of the 
time-dependent inviscid flow properties will be required. 
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4. The Flow in the Viscous Layer 


4.1 Governing Equations 


The viscous region consists of thin boundary layers that lie along the upper and lower 
surfaces of each blade and thin wakes which extend downstream from the blade trailing 
edges, as indicated in Figure 3. We assume the blade motions to be of sufficiently small 
amplitude, i.e., on the order of the viscous-layer displacement thickness (£), that the system 
of field equations that governs the flow in a viscous layer is invariant under a transformation 
from space- to blade-fixed coordinates. Thus, in terms of the blade-fixed r, n coordinates we 
can write 


Sp d(pU) djpV) 

dt dr dn 



(4.1) 

(4.2) 


and 


~DH dP wl d 
P Dt ~ dt + ^ dn 


(l - (TV) -1 ) U ~ + 


dH 


dn Pr dn 


— ph'v' 


(4.3) 


Here, U and V are the velocity components, relative to the blade-fixed frame, in the positive 
t- and n-directions, respectively, and p is the molecular viscosity of the fluid. The terms 
that account for the acceleration of the blade-fixed frame are assumed to be small and have 
therefore been neglected in reducing the Navier-Stokes equations to viscous-layer equations. 

In equations (4.1)-(4.3) the coordinates r and n measure distances along and normal to, 
respectively, a blade surface or a reference wake streamline, and D j Dt = dfdt -f V • V is 
the convective derivative operator in the blade-fixed frame. The reference wake streamline 
emanates from a blade trailing edge and is entirely contained within the viscous wake. A 
distinction should be made between the independent variables that describe the flows in the 
upper and lower surface boundary layers and in the wakes, e.g., by attaching subscripts to 
r and n. However, as a convenience, we are neglecting to make this distinction explicitly. 

The symbols p , U, V, P, H, and p in Eqs. (4.1) — (4.3) refer to ensemble (or Reynolds) 
averaged values of the fluid dynamic variables; u\ v' and h' are the values associated with 
random turbulent fluctuations, and the overbar indicates a turbulent correlation, which 
must be determined empirically. As a consequence of the high Re and small-amplitude 
approximations, the pressure in the thin viscous layer is a function only of r and t, and 
H = T + U 2 / 2. Since the interaction between the flows in the viscous layer and the external 
inviscid stream is assumed to be weak, the pressure and the flow properties at the edge of 
the viscous layer are determined by the inviscid values of these variables at the blade and 
wake surfaces. 

In addition to the foregoing conservation equations, the equation of state for a perfect 
gas, i.e., 

P = , (4.4) 
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a law relating the molecular viscosity to the temperature, e.g., 



T \ 3/2 7Lqq ±Tc 
T- oo) T + T c ’ 


(4.5) 


and equations relating the turbulent correlations u'v' and h'v' to the ensemble- averaged flow 
quantities, i.e., 

dU . dH t dH 


pu'v^-e--^ and ph'V = ^ , 


(4.6) 


are also required. Equation (4.5) is a form of Sutherland’s Law. Here ji-oo is the molecular 
viscosity at the temperature T—oo &nd Tc is 3- constant, which for air has a dimensional 
value of 110°K [24]. The turbulent correlations u'v 1 and h'v' have been related to mean flow 
gradients, using Prandtl’s mixing-length hypothesis. Here e and tfj are the eddy viscosity 
and eddy diffusivity, respectively, and play roles similar to their molecular counterparts. 
The eddy diffusivity is related to the molecular diffusivity through the introduction of the 
turbulent Prandtl number, Prj- 


Initial and Boundary Conditions 

The foregoing system of field equations is parabolic in time and in the streamwise direc- 
tion. Therefore, the streamwise component of the velocity and the total enthalpy must be 
known for all time at some upstream streamwise location, and these variables, along with 
the normal velocity, must be known throughout the solution domain at some initial time. 
Also, conditions on the fluid properties at the edge(s) of the viscous layer, i.e., 

U->U e {T,t) and H —> H e (r,t) for n — ► Too , (4.7) 

where the limits Too and — oo refer to the edges of the upper (+) and lower (— ) surface 
boundary layers and the upper and lower edges of a wake, a no-slip condition and either a 
prescribed temperature or heat flux condition at a solid blade surface, i.e., 

dH 

U = v = 0 and H = H w (t, t) or — = Q w (r,t) for n = 0 , r < t u ; (4.8) 

and a condition on the fluid velocity normal to a reference wake streamline, i.e., 

V = 0 for n = 0 , t > Tie , (4.9) 

must be enforced. Here the subscripts w and e denote the values of the fluid properties at a 
solid wall and at the edge of the viscous layer, respectively, and the subscript te refers to the 
airfoil trailing-edge point. The fluid velocity, U e , and total enthalpy, H e , at the edges of the 
viscous layers are determined by the inviscid solution along the blade surfaces and reference 
wake streamlines [c.f. Eqs. (3.44) and (3.48)]. 


21 


Turbulence and Transition Models 


The models used here and in Ref. [13] to simulate the effects of turbulence and transition 
on the flow in the viscous layer are the algebraic eddy- viscosity model proposed by Cebeci and 
Smith [25], the transition length correlation proposed by Dhawan and Narashima [26], and 
the wake turbulence model proposed by Chang et al [27]. Also, since flows in turbomachines 
are known to be characterized by high freestream turbulence levels, a simple modification 
developed by Yulias [28] has been incorporated into the turbulence model to account for the 
effects of freestream turbulence on the viscous layer. These models are easy to implement, 
and are known to be reasonably accurate for steady flows with mild pressure gradients. Their 
ability to accurately represent turbulence and transition in unsteady flows is not known, 
however; therefore, results must be interpreted with some caution. 

The Cebeci-Smith algebraic model divides a solid-surface boundary layer into inner and 
outer regions, where e = c, and e = e D , respectively. The inner model is applied from the 
wall out to the point at which e, = e 0 ; the outer model, from this point to the edge of the 
boundary layer. The eddy viscosity in the inner region is given by 


e, = 7 t r p( 0 . 41 n) 2 [l — exp(— n/Ax)] 2 Re 


dU_ 

dn 


where 


At = 


26/z 
pV Re 




P P 


dn 


1 - 11.8P+ 


fpw\ fhV ) 

\h) Khz J J 


P + = v/R^ 

p 


k ( du . 


-i 


dt 


+ 


fT dUA (~-x--dUY 2 

Ue dT)[ P fi dn) w ’ 


and t tT is a longitudinal intermittency factor which models transitional flow. 
The eddy viscosity in the outer region is given by 


r°° ~ ~ 

Eo = ItrRepUeX / (1 - U/U e )dn , 

Jo 


(4.10) 


(4.11) 


(4.12) 


where 


and 


X = 1.55xo(l + 7r) 1 , 7r = 0.55[1 — exp(— .243zl^ 2 — .298zi)] , 
f Re#/425 — 1 for Re tf > 425 
z i = \ 

{ 0 for Re# < 425 , 


Re^ 


pJjjte r 
K Jo u e { 



(4.13) 


The parameter tt heis been defined to account for low momentum thickness Reynolds number. 
The Clauser constant, Xo, is usually set equal to 0.0168, but following Yuhas [28] we set 


Xo = 0.0168[1.0 + 18.4(v' • v'^ + 99.6(v' • v')-oo] , 


(4.14) 
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where (v' is the mean-square of the freestream turbulence velocity, to account for 

the effects of freestream turbulence on the development of the turbulent boundary layer. 

Blade boundary layers can contain transitional regions of significant extent. Therefore, 
a transition model is needed to properly predict boundary-layer flows. Several models have 
been proposed, based on correlations with experimental data, for steady transitional flows. 
One of these is the Dhawan-Narashima forced transition model [26] in which the intermit- 
tency factor is given by 


7tr = 1 ~ exp 



(4.15) 


Here and r 2 are the streamwise locations at the beginning and end of the transition region, 
respectively. These locations can be specified as functions of time to model unsteady flows. 

The eddy viscosity in the wake is based on the model developed by Chang et al [27]. At 
each streamwise station in the wake Equation (4.13) is evaluated twice: once for — oo < n < 
n min> where is the location at which the streamwise velocity reaches a minimum value, 
and once for n^n < n < oo. The maximum of these two values, e max , is then used to set the 
eddy viscosity, eyv, at each streamwise position in the wake, i.e., 


CyV — fnui [fnux fmux.te] exp 


T ~ T te 


206, 


tc -» 


(4.16) 


In this equation r — r te is the distance measured along the reference wake streamline from 
the trailing edge point, and 8^. is the sum of the upper and lower surface boundary-layer 
thicknesses at the trailing edge. The boundary- layer thickness is defined to be the normal 
distance from a blade surface to the point at which U /U e = 0.995. 


4.2 Transformed Equations 

For laminar flows, the Levy-Lees transformation (see Blottner [29]) provides independent 
variables that effectively capture the growth of the viscous layer with increasing streamwise 
distance. In addition, the transformed equations reduce to similarity equations at an airfoil 
leading edge. These features facilitate the determination of a numerical solution. Thus, an 
extension of the Levy-Lees transformation, in which the laminar edge viscosity is replaced 
by an effective turbulent viscosity [30], is applied herein to accommodate unsteady turbulent 
flows. 

We define new independent variables £ and i] according to 

£ = J q Qdr and tj = j Q P dn » ( 4 - 17 ) 

where Q = p e U e p e [ 1 + (e 0 / p) e ] and l and /2 are the turbulent eddy and molecular viscosities 
in the steady background flow. In Ref. [13] the temporal mean values of the fluid properties, 
p e , U e , and p e , at the edge of the viscous layer were used to define Q ; however, in the present 
effort in which the unsteady inviscid flow is regarded as a small perturbation of a nonlinear 
mean flow, we can replace the temporal mean values of these variables by their values in the 


23 


steady background flow to avoid introducing additional nomenclature. The coordinate £ in 
(4.17) is a function of r alone; the coordinate T], by contrast, is a function of r, n, and t. 
We also introduce the new dependent variables 



n 

Q 


\U e dt dr 



and Q 



(4.18) 


where Q = pfi&A l + (e/^) e ]. After substituting Eqs. (4.17) and (4.18) into the viscous-layer 
field equations, (4.1), (4.2) and (4.3), and performing some algebra, we find that 


2 ^^ + |y + a _ 3 = 0> 

Q S( dn Q 


(4.19) 
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where the following parameters have been introduced: 


(4.21) 
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Pehi 1 + (t/p)e) ’ 


(4.22) 


Special attention must be exercised in applying the foregoing transformation to unsteady 
wakes. In the present analysis all quantities on the upper (i.e., n > 0) and lower (n < 0) sides 
of a wake are referenced to their respective edge conditions, i.e., we set^^F = U jU+ for n > 0 
and T = U /U~ for n < 0. In general, the inviscid streamwise velocity, U e , and total enthalpy, 
H e , w iH be discontinuous across unsteady wakes. As a consequence, the dependent variables 
T and Q are discontinuous across reference (n = 0) wake streamlines. These discontinuities 
must be taken into account in developing a numerical solution procedure, to ensure that the 
physical variables (U, V and H) are continuous across viscous wakes. 
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Also, the upper- and lower-surface values of the independent variable £ will generally 
differ at a blade trailing edge. Therefore, in the present effort the upper-surface ^-variable 
is used to continue the viscous calculation into the wake, i.e., 

6v = £ + r Q + dr w . (4.23) 

Jrt* 

However, two different rj- coordinates, i.e., 

Vw = j Q ~P dn ' n <°> (4-24) 

are used for the wake calculation. Thus, at the trailing edge, the lower-surface boundary- 
layer solution must be interpolated onto the wake rj mesh before continuing the viscous 
calculation into the wake. 

Boundary and Initial Conditions on T and Q 

In addition to satisfying the field equations (4.19) — (4.21), the dependent variables T 
and Q must be prescribed functions of £ and i] at some initial time, and prescribed functions 
of tj and r at initial or upstream streamwise locations on the upper and lower surfaces of 
each blade. Also, in terms of the transformed variables, the following conditions must be 
imposed at the edge(s) of the viscous layer, on the blade surface, and along the reference 
wake streamline, c.f. Eqs. (4.7), (4.8) and (4.9): 



T — ► 1 and Q — > 1 for |r/| — > oo , 

(4.25) 


T = 0, V = 0, and Q — Q w (£, r) or 



t nr* / o \ “(2£/He)a 4 L(£,t) r n 

( dg/dTj ) w = - - for T) = 0, ( < , 

P e U e H e 

(4.26) 

and 

V = 0 for r) = 0, ( > , 

(4.27) 

respectively. 




For the calculations reported herein, the condition at the initial time is the viscous-layer 
solution in the absence of unsteady excitation, i.e., the solution for steady background flow. 
The calculated transients resulting from this approximation have been found to die out with 
increasing time. The upstream profiles, required at each time step of the unsteady viscous- 
layer calculation, are obtained from a similarity analysis of the flow in the stagnation region, 
as described below. 

4.3 Stagnation Region 

We seek a similarity solution for the flow in the stagnation region, which can be used 
to set the instantaneous upstream profiles, and required to initiate 

a full viscous-layer calculation. The overall strategy is to develop a similarity solution for 
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a simple flow configuration and then, to adapt this result for application to the unsteady 
cascade problem. The present analysis for unsteady, compressible flow is modeled after 
the incompressible analyses developed independently by Rott [31] and Glauert [32]. Their 
analyses provided an exact solution to the Navier-Stokes equations for incompressible flow. 
This is not the case for compressible flow, however, because it is necessary to neglect the 
dissipation term in the energy equation to obtain a similarity solution. Fortunately, frictional 
heat dissipation is of limited importance for a wide range of practical flows. 

We consider two-dimensional compressible flow around a flat plate, which is oriented 
normal to the stream direction and undergoes a harmonic motion at velocity iur T exp(i'u>f) 
in its own plane. As a convenience to this discussion, we use a complex-variable description 
for the unsteady flow properties. It is to be understood, however, that the real parts of 
the various complex parameters represent the actual physical variables. The flow in the 
viscous layer will be determined in terms of spaced-fixed Cartesian coordinates, Tg and 
jig. The resulting stagnation-region solution can then be expressed in terms of plate-fixed 
coordinates as a prerequisite to its implementation into the complete unsteady viscous-layer 
calculation. The Tg and ng axes are directed along and normal to the plate, respectively, 
and the coordinate origin lies at the point at which the dividing inviscid streamline impinges 
on the plate. The fluid velocity components in the Ts- and ng-directions are Us and Vs, 
respectively. 

The inviscid flow is steady relative to the space-fixed frame, with velocity and pressure 
gradient at the edge of the viscous layer (i.e., at n s = 0 in the inviscid region) given by 

V S ,e = Us,e e T S = CTge^ , (4.28) 


and 


VP = 



-KUs„ 


dU 


SjC 


dr* 


~T S 


p c ^ ^se rs 


(4.29) 


The constant c in these equations is determined by the behavior of the inviscid solution in 
the vicinity of the mean stagnation-point location (r s = 0). The total enthalpy at the edge 
of the viscous layer is constant, i.e., H e = T e + f7| e / 2 = //_<»• 

The flow in the viscous layer is described by the field equations (4.1)— (4.5), the edge 
conditions (4.7) and the surface conditions (4.8), except that in the present case 


Us(rs,ns,t) = iu>r T exp(iwt) for n s = 0 . 


(4.30) 


We assume that this flow is laminar and that the heat generated by viscous dissipation is 
negligible, i.e., ftd 2 U s /dnl « 0. The streamwise velocity and temperature can be expressed 
in the forms 

U s (t s,V,i) = cTsf'iv) + ivr T g(r])e ,wt , (4-31) 


and 

T{rs,Ti,t) = f e (T S ,t)0(r}) , 


where 



(4.32) 

(4.33) 
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Note that dP/dns ~ 0 and Vs <C Us] therefore, p e /~p = T jf e — 9 and H = T + U$/2 in the 
viscous layer. 

After substituting the foregoing results into the field equations that govern the flow in the 
viscous layer and performing the necessary algebra, we obtain an expression for the normal 
component of the fluid velocity, Vs, and the following set of ordinary differential equations 
for /, g and 9. 

f' 2 -fr-(lf")'-0 = 0 , (4.34) 

(- + f)g - fg' - W = o , (4-35) 

c 

and , 

(^) + f0, = 0 ’ ( 4 - 36 ) 

where / = ppf(p e fi e ). The conditions on /, g and 0 at the plate surface (g = 0) and at the 
edge of the viscous layer (77 — > 00 ) are 

/( 0) = /'( 0) = 0 , 5 (0) = 1 and 9'( 0) = 0 or 0(0) = T w /T e ; (4.37) 

and 

/'( 00 ) = 1 , 5 ( 00 ) = 0 and 0(oo) — 1 . (4.38) 

Note that for steady (g = 0), incompressible (0 = 1) flow, we recover the classical stagnation- 
region problem studied by Hiemenz. Also, Eq. (4.36) is identical in form to the energy 
equation that describes steady compressible stagnation-point flow [33]. Finally, the present 
analysis can be readily extended to consider more general unsteady motions, e.g., 

N 

Us = CTsf'iv) + 5Z iu nr T ,ngn(g)exp(iu} n t) . (4.39) 

n=l 

In this case solutions for N uncoupled equations for the g n , n = 1, 2, 3, . . . , A r , are required. 

The set of nonlinear ordinary differential equations for /, 0 and g can be solved using 
an implicit finite-difference technique similar to that used to obtain nonsimilar solutions to 
the full unsteady viscous equations (cf. §4.4). The functions / and 0 must be determined 
simultaneously, but they are independent of g , which can be determined once / is known. 

To provide the upstream profiles required for a complete viscous-layer calculation it is 
convenient to express the stagnation-region solution in terms of body- or plate-fixed coordi- 
nates. Thus, we consider Cartesian coordinates t, n, where the r- and n-axes are parallel to 
the spaced-fixed t$- and ns-axes, respectively, and the mean position of the r, n coordinate 
origin coincides with the origin of the spaced-fixed Tg, rig-frame. The streamwise positions 
and velocities of a fluid particle in the two coordinate frames are related by 


rg = r + r T exp (iu>t) 

(4.40) 

and 


Us = U + iur T exp(iut) . 

(4.41) 

After combining (4.31), (4.40) and (4.41), we find that 


U = erf '(g) + b\(g) exp(iut) , 

(4.42) 
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where 


b = cr T (l — iuj/c) , 


and 

(1 - iu/c)\(rj) = /'(f?) + (iu}/c)\g(rj) - 1] . 

At the plate surface U — 0, and at the edge of the viscous layer 

U(rj — > oo) = U e = ct + bexp(iu>t) . 

The total enthalpy H = t/|/2 + T is given by 

H = (ct/') 2 /2 + c( 1 — ^a;/c) -1 6r/ , (/ , + iu>g/c)exp(iojt ) 

+(1 — iu/c)~ 2 b 2 (f + iug/c) 2 exp(2iu;t)/2 + T e 0 . 


(4.43) 

(4.44) 

(4.45) 


(4.46) 


The locations, i.e., r = ±|tj|, at which the stagnation-region solution can be applied as 
an initial condition for an unsteady viscous-layer calculation are determined by two criteria. 
First, the stagnation-point motion must be contained within the interval [— |t/|, |r/|] and 
second, to avoid stability problems in the subsequent boundary-layer calculation, those points 
at which flow reversals, i.e., ( dU / dn) w < 0, occur must also be contained within this interval. 
Reverse flow is associated with the lag in the response of the low-momentum fluid near the 
wall to changes in the velocity at the edge of the viscous layer. The extent of the interval 
over which the stagnation point moves, i.e., [— \tm\, \tm |], is determined from the inviscid 
velocity distribution (4.45) and is given by [— c _1 |6|, c|6|]. The extent of the reverse-flow 
interval, i.e., [— |tr|, |tr|], is determined by the maximum and minimum values of r for 
which (dU/dr}) w = 0. We find that 


I t h| < c 1 \b\ 


i + wc-y(o)//"(o) 

1 — iui/c 


(4.47) 


Thus, |r/| > max(|rjvf|, |th|). For our application to unsteady cascade flows, we have found 
that setting |r/| = 2c -1 |6| leads to reasonable results in the stagnation region and allows us 
to continue the viscous-layer calculation over blades and their wakes. 

In applying the foregoing stagnation-region analysis to unsteady cascade flows, we use 
the inviscid velocity distribution along a blade surface, i.e., 


t) = U c (t ) + UeO") exp(iuf) -1- . . . , 


(4.48) 


where U e (r ) and u e (r) are the nonlinear steady and the complex amplitude of the linearized 
unsteady velocities, respectively, to determine the parameters c and b in (4.42) and (4.46). In 
particular, these parameters are determined by matching the analytical velocity profile (4.42) 
to the calculated inviscid velocities at ±|r/|, i.e., we set 


c = (<9£/ e /dr) T= o and b = u e (±r/) 


(4.49) 


Thus, with u> specified as an input quantity, the functions / 0 and g determined by an 
implicit finite-difference procedure and the parameters b and c given by (4.49), the streamwise 
velocity (4.42) and total enthalpy (4.46) profiles needed to initiate the unsteady viscous-layer 
calculation can be specified. 
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4.4 Numerical Solution Procedure 

The transformed form of the viscous-layer equations [Eqs. (4.19)-(4.21) and (4.25)-(4.27)j 
are solved by using a finite- difference approach to approximate the various partial deriva- 
tives and to convert the system of differential equations into a system of algebraic equations. 
Streamwise and temporal derivatives are approximated by first-order backward differences 
and normal derivatives are approximated by second-order central differences. The stream- 
wise momentum and energy equations contain nonlinear terms, and therefore, they must be 
linearized to facilitate their solution at each point in space and time. This is accomplished 
by using a Newton iteration in which the initial guess for the profiles comes from the solution 
at the previous spatial or temporal step. The equations are solved using local iteration to re- 
move the linearization error, and repeating the iteration until the values of the flow variables 
converge to within a specified tolerance level. The finite-difference approximation results in 
a block tridiagonal system of linear algebraic equations at each step, which is solved using 
the Thomas block inversion algorithm. 

The analyses for surface boundary layers and wakes are almost identical, except for the 
implementation of different boundary conditions. In addition to replacing the surface no- 
slip condition with a zero normal velocity condition at the reference wake streamline, the 
possibility of jumps in the tangential velocity and total enthalpy across a viscous wake must 
be taken into account. As discussed previously, the transformed governing equations are 
written in terms of two different sets of variables, each applying on one side of the wake 
streamline, and scaled to the edge conditions for that side. For example, on the upper side 
of the wake, F + = u/u+ and G + = H/H+ ; on the lower side, F~ = u ju~ and G~ - H/ H~ . 
The discretized equations on either side of the reference wake streamline are written in terms 
of that side’s variables. However, at the wake streamline, (rj = 0), variables from both sides of 
the wake are used in the momentum and energy equations, due to the use of central-difference 
approximations for the ^-derivatives. For consistency, the equations must be written in terms 
of a single set of variables. This is accomplished here by writing the equations in terms of the 
upper surface quantities. Thus, whenever a lower surface variable appears in the equations 
at r) = 0, it is rewritten in terms of the upper surface edge conditions. For example, the 
variable F~ is written in upper variables as F~ • (u~ /u+) = u - /uj". Further details on the 
reference wake streamline procedure are given in the Appendix, where the finite-difference 
solution procedure for the viscous layer is discussed in detail. 

The finite-difference approximation is implemented on a grid that is nonuniform in both 
spatial directions, and uniform in time. The grid distribution in the streamwise direction 
is chosen so as to cluster points near the blade leading and trailing edges, with each mesh 
interval being set equal to a constant times the previous one, according to the relation 

A6+i = ^A6, / = 1, 2, ... ,X (4.50) 

where l is the streamwise mesh point index (/ = 1 at the blade leading edge and increases 
with distance downstream), > 1 from the leading edge of the blade to midchord, and 
< 1 from midchord to the trailing edge. The grid is then stretched (i.e., > 1) aft of 

the trailing edge, and generally extends between one and two blade chords downstream of 
the trailing edge. 
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A stretched grid is also employed in the normal direction, with clustering near the surface 
to capture the large velocity gradients that occur near the wall, and to ensure that, for 
turbulent flow, there are enough points in the near-wall region to adequately resolve the 
laminar sublayer. Geometric stretching is again employed, with 

At7„ +1 = K v AT) n , n = l, 2, ... ,N (4.51) 

where n is a normal mesh point index (n = 1 at the blade surface or reference wake streamline 
and increases with increasing distance into the flow). Since the resolution requirements are 
a function of the solution, the normal grid stretching, K v , and the spacing at the wall, Ar}i, 
must be determined for each case. For turbulent flow, this is accomplished by monitoring 
the value of Y + at the first mesh point from the wall (i.e., F^t. 2 ), where Y + is a Reynolds 
number based on the friction velocity, {t w Ip w ) 1, ‘ 2 and the normal distance from the blade 
surface. For accuracy, should be close to one, which can be achieved by adjusting 
The stretching parameter K n is then adjusted to place the outer edge of the mesh, 77 = 77^, 
far enough from the surface to allow the edge conditions to be approached asymptotically. 
The number of points used in the normal direction is chosen to allow the flow over the entire 
viscous layer to be resolved accurately at all streamwise stations — generally, between 50 
and 100 points are sufficient. 

The constant value of the time-step used for the temporal discretization is chosen based 
on the frequency of the imposed unsteady excitation and the number of time-steps prescribed 
for each period of oscillation, generally between 20 and 50. 
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5. Numerical Results 


Response predictions will be given below for two-dimensional blade rows operating at 
subsonic Mach numbers and subjected to external unsteady aerodynamic excitations to 
demonstrate the present capabilities of the linearized inviscid unsteady flow analysis, LIN- 
FLO, and the unsteady viscous-layer analysis, UNSVIS. The inviscid analysis will be applied 
to predict the response of compressor and turbine cascades, and their flat-plate counterparts, 
to vortical excitations; the unsteady viscous analysis, to study the viscous-layer responses of 
flat-plate and turbine cascades to incident acoustic excitations. 


5.1 Inviscid Flow: Vortical Excitation 


We consider two of the cascades studied in Ref. [16] — a compressor exit guide vane 
(EGV) and a turbine cascade proposed as the fourth standard cascade configuration for 
turbomachinery aeroelastic investigations by Fransson and Suter in Ref. [34]. In addition, 
we will compare the unsteady responses of these compressor and turbine configurations to 
those of corresponding flat-plate cascades, i.e., flat-plate cascades having the same blade 
spacing, G , and operating at the same inlet flow conditions, M_oo and fl-oo, as the real- 
blade configuration. However, the flat-plate blades are staggered so that the blade mean 
positions are aligned with the mean inlet flow direction, i.e., 0 = Therefore, the local 

steady Mach number, M = and flow angle, fi = are constants throughout the 

flat-plate flow fields. The steady background flows through the EGV and turbine cascades 
have been determined using the methods of Ref. [18]. In each case a Kutta condition has 
been applied at blade trailing edges; therefore, only inlet uniform flow information, e.g., 
M_ <X) and fi_oo> must be specified for the steady calculations. 

First-harmonic unsteady solutions have been determined, using LINFLO, on a global H- 
type mesh and, for one EGV calculation, on a composite mesh constructed from this H - mesh 
and a polar-type local mesh. The H - mesh used here consists of 40 mean-flow streamlines 
and 155 “axial” lines and extends one axial chord upstream and downstream of the blade 
row. The mesh streamlines are packed near the blade and wake surfaces; the axial lines, near 
the blade edges. The local mesh used in the composite-mesh EGV calculation consisted of 
100 radial lines and 11 “circumferential” lines, with the latter being tightly packed near the 
blade surface. As part of this study, classical linearized theory (CLT) response predictions 
have also been determined for the flat-plate cascades using the analysis of Smith (Ref. [35]). 

The unsteady flows are excited by vortical gusts. The temporal frequency, u>, interblade 
phase angle, a, and complex amplitude of the gust velocity component normal to the 
freestream direction, i.e., Vr- 0 0 • e^v, are prescribed quantities. The component of the gust 
velocity in the inlet freestream direction, Vr _<» • er, is determined by the divergence-free 
condition, i.e., it • vr-oc = 0. The wave number #c_ has a component = aG~ l 

in the cascade- or ^-direction and a component k t = — in the inlet freestream 

direction. Therefore, 


K_ 


oo 


— (u> sec fi-oo + vG 1 tan + oG 

— ujej + (wtanfi-oo + oG~ x secf)_oo)ejv . 


(5.1) 
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(5.2) 


For a pure vortical excitation the rotational velocity field is given by 

v H (x) = V(X • v R ,_oo) exp(i/c_ 00 • X) , 

where X is defined in Eq. (3.18). For the flat-plate cascades and far upstream of the EGV 
and turbine blade rows v/?(x) = ~vr,-cx> exp(iK_oo • x), and the complex amplitudes of the 
gust vorticity is C-oo = * K -<» x V R- oo- 

Since the inlet and exit flows considered here are subsonic, a subresonant unsteady ex- 
citation produces an acoustic response in which all waves attenuate with increasing axial 
distance from the blade row. A superresonant (m,n) excitation produces a response in 
which m and n such waves persist far upstream and/or far downstream, respectively, and 
carry energy away from the blade row. An acoustic resonance occurs if at least one acoustic 
response wave persists in the far field and carries energy along the blade row. 

We will examine the response of the compressor, turbine and corresponding flat-plate 
cascades to vortical gusts with Vr,-oo • e# = (1 5 0), u> = 5 and — 37r < a < — 7r. It should 
be noted that u> and a must be of opposite sign to model a realistic cascade/ vortical gust 
interaction. In the present study we set u> > 0, and therefore, a must be less than zero. 
Theoretical predictions will be given for the Mach number, drift function and stream function 
in the steady background flows, and for the vorticity, source term, velocity potential and 
pressure in the linearized unsteady flows through the compressor and turbine cascades. We 
will also examine the behavior of the unsteady pressure differences [Ap(x) = p(xg) - p(xj)] 
and unsteady aerodynamic moments (m) acting on the reference (m = 0) blades of the 
compressor, turbine and corresponding flat-plate blade rows. 

Compressor Exit Guide Vane 

The blades of the compressor exit guide vane (EGV) are constructed by superimposing 
the thickness distribution of a NACA 0012 airfoil on a circular arc camber line with height 
at midchord of 0.13. This cascade has a stagger angle 0 of 15 deg, a blade spacing G of 0.6 
and operates at an inlet Mach number and inlet flow angle of 0.3 and 40 deg, respectively. 
The calculated exit Mach number, exit flow angle and mean lift force, F y e y , acting on each 
blade are 0.226, -7.4 deg and 0.36, respectively. The predicted steady Mach number field 
and Mach number distribution along a blade surface for the steady flow at M-<x> = 0.3 
and fl_oo = 40 deg through the EGV are shown in Figure 4; the drift and stream function 
contours, in Figure 5. Far upstream of the blade row, the drift function contours are parallel 
to each other and perpendicular to the inlet freestream direction, but as the flow proceeds 
downstream, gradients in the steady velocity produce distorted contours. These distortions 
are mild over most of the blade passage, but quite severe in the immediate vicinity of the 
mean blade and wake surfaces, where the drift function contours stretch downstream from 
the leading-edge stagnation point. 

The meshes used to determine global- and composite-mesh solutions for the unsteady po- 
tential are depicted in Figure 6. The streamline global mesh shown here was also used to eval- 
uate the rotational velocities, \r and v,, the vorticity, £, and the potential-equation source 
term, p _1 V • (pv*), over the extended blade passage solution domain. For the composite- 
mesh potential calculations, the values of these quantities at the local-mesh points were 
determined by an interpolation from the global to the local mesh. 
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Global- and composite-mesh solutions for the in-phase component (real part) of the 
unsteady potential due to a vortical excitation at u = — 2x are shown in Figure 7 the global- 
and composite- mesh results for the unsteady pressure-difference response, in Figure 8. The 
results in Figures 7 and 8 indicate that there are only slight differences between the global 
and composite-mesh unsteady flow solutions, and hence, that the streamline H- mesh of 
Ref. [12] and [16] is adequate for predicting low-speed unsteady cascade flows. The global 
mesh analysis is very efficient; for example, a complete unsteady flow calculation on a 155 x 40 
global mesh requires about 2 1/2 minutes of CPU time on an APOLLO 10000 work station. 
Most of the information determined from such a calculation can be saved and re-used in 
subsequent unsteady flow calculations for different frequencies, interblade phase angles, and 
modes of excitation, which then require only about 23 seconds of CPU time for each unsteady 
case. Composite- mesh calculations typically require four (4) times the CPU time required for 
a global mesh calculation. It should be noted, however, that the composite-mesh analysis 
represents an important solution capability because it allows a great deal of flexibility in 
prescribing a mesh structure for accurately analyzing linearized unsteady flows. In addition, 
this type of capability will be essential for analyzing transonic problems in which a mesh 
that conforms to the mean shock loci is required to accurately predict the effects of shock 
displacements. 

We proceed to examine in somewhat more detail the unsteady flows through the EGV 
that are excited by vortical gusts with vr-oo'^n = (1,0) and w = 5. These results (Figures 9 
through 15) and similar ones for the turbine in Figures 18 through 24 were determined using 
global mesh calculations. Contours of the in-phase component or real part of the unsteady 
vorticity, source term and pressure are shown in Figures 9, 10 and 11, respectively, for vortical 
excitations at a = — 7r, — 2 t and — 3 tt. The wave-number magnitudes, |k_oo|, associated 
with these gusts are 5.66, 10.73 and 17.08, respectively, and the arguments relative to the 
axial flow direction, a = tan _1 (/c, i _ 00 /«: {( _ 00 ), are -112.1 deg, -77.8 deg and -67.0 deg, 
respectively. The vortical gusts are distorted as they are convected by the nonuniform mean 
flow through the EGV. The vorticity contours in Figure 9 and the source term contours in 
Figure 10 indicate that this distortion increases in severity, i.e., the vorticity and source-term 
contours are more severely stretched and re-oriented within the blade and wake passages with 
increasing values of |k._oo|. The results in Figure 10 reveal the rather strong variations in the 
source term that occur throughout the extended blade passage solution domain, particularly 
for the high wave number gust. 

The vortical gusts at a = —it and a = — 3zr produce a subresonant pressure response, 
i.e., all acoustic response waves attenuate with increasing axial distance from the blade row 
(see Figure 11). The vortical excitation at <r = —2 7r produces a superresonant (1,1) response 
in which a pressure response disturbance at lo = 5 and k v = 0 propagates away from 
the blade row in both the far upstream and far downstream directions. Note the spurious 
pressure behavior indicated in Figure 11 at the upstream boundary (for <r = —2 tt) and at the 
downstream boundary (for a = —27 r, — 37r) of the solution domain. This can be attributed 
to the nonuniform spacings between the mesh streamlines currently used in the global mesh 
calculations. The global mesh streamlines must be packed near the blade and wake surfaces to 
resolve the vortical flow near these surfaces, but this leads to a reduced accuracy in matching 
the numerical near-field and analytical far-field solutions for the unsteady potential at the 
inlet and exit boundaries, and thus, in predicting the acoustic response in the far field. The 
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use of a composite- mesh procedure in which the global mesh contained uniformly spaced 
streamlines might overcome this problem. 

Responses at the reference blades of EGV and corresponding flat-plate cascades are 
depicted in Figures 12-15. The flat-plate cascade has a stagger angle of 0 = ft = 40 
deg, a blade spacing G = 0.6, and operates at a Mach number, M, of 0.3. The flat-plate 
predictions presented in Figures 12-15 have been determined using both LINFLO and the 
classical linearized analysis of Smith (Ref. [35]). Recall that in the general linearization 
the unsteady potential equation (3.11) contains a source term, i.e., p -1 V • (pv*), and the 
normal derivative of the potential, • n, is zero at the mean blade surfaces. For flat- 
plate blades aligned with the inlet flow direction the source term reduces to V 2 </>*. In the 
classical linearization the potential equation is homogeneous, and the normal derivative of 
the potential cancels the normal component of the gust velocity at the mean blade surfaces. 
Thus, a comparison between the present and the classical linearizations is quite meaningful 
because it reveals the ability of the numerical field methods, associated with the former, to 
accurately account for the strong source term variations that occur over the solution domain. 

The pressure difference distributions that act on the reference blade of the EGV and flat- 
plate cascades for vortical excitations at V/j,_ «, • — (1,0), u> — 5 and cr = — x, cr = —2 x 

and a = — 3x are shown in Figures 12, 13 and 14, respectively. The unsteady moments 
acting on the reference blades of the two cascades are plotted versus interblade phase angle 
for vortical excitations at Vr-oo'^n = (TO), w = 5 and — 3x(— 540 deg) < a < — x(— 180 deg 
in Figure 15. Here, the excitations at a — —404.2 deg and —293.9 deg produce resonant 
acoustic response disturbances far upstream and far downstream of the flat-plate cascade and 
far upstream of the EGV; those at cr = —414.3 deg and —308.8 deg produce such responses 
far downstream of the EGV. 

Vortical gusts are convected without distortion by the uniform mean flow through the 
flat-plate cascade, i.e., the rectilinear vorticity contour patterns that exist far upstream of 
the EGV blade row are maintained throughout the flat-plate flow field. The results in Fig- 
ures 12-14 indicate the effects of gust distortion, due to nonuniform mean flow phenomena, 
on blade pressure- difference responses. For the vortical gusts at cr = — x and a = — 2x the 
pressure-difference distributions along the EGV and flat-plate blades show somewhat similar 
qualitative behaviors, but important quantitative differences. Such differences between real 
blade and flat-plate responses coincide with intuitive expectations. In contrast, the responses 
of the EGV and flat-plate blades to the vortical gust at a = — 3x bear almost no qualita- 
tive resemblance. This result is somewhat surprising and, if correct, indicates the effect 
of mean-flow nonuniformity on the unsteady pressure response of cascades to high-wave- 
number vortical excitations to be quite substantial. The real blade and flat-plate moment 
responses shown in Figure 15 again show some qualitative similarity except, perhaps, over 
the phase-angle range —540 deg < cr < —420 deg. 

Note that the LINFLO and Smith analysis (CLT) pressure-difference predictions in Fig- 
ures 12 through 14 for the flat-plate cascades are in excellent agreement. Indeed, the pressure- 
difference curves predicted by the two analyses are almost coincident for the gusts at cr = — x 
and a = — 2x, and only slight differences exist for cr = —3 x. The LINFLO and classical the- 
ory predictions for the moment response (Figure 15) are in excellent agreement over the 
entire interblade phase angle range, i.e., -3x < cr < -x. 


34 



Turbine Cascade 


As a representative turbine configuration we have selected the fourth standard config- 
uration of Ref. [34], but for the present study we have modified the blade profiles defined 
in Ref. [34] so that our example blades close in sharp, i.e., wedge-shaped, trailing edges. 
The turbine cascade and the corresponding flat- plate configuration operate at an inlet Mach 
number of 0.19 and an inlet flow angle of 45 deg, and the blade spacing G is 0.76. The 
turbine cascade has a stagger angle of 56.6 deg, and the flat-plate blades are staggered at 
0 = O = 45 deg. The calculated exit Mach number and flow angle for the turbine are 0.49 
and 72.0 deg, respectively, and the mean lift acting on each blade is —2.11. The predicted 
steady Mach number contours and blade-surface Mach number distributions for the steady 
flow at Af.oo = 0.19 and f !_<*, = 45 deg through the turbine are shown in Figure 16. The 
drift and stream function contours for this flow are shown in Figure 17. Note the stretching 
of the drift function contours within the blade passage, which is caused by the acceleration 
of the steady flow. 

Linearized unsteady flow predictions for the turbine cascade, as determined using the 
LINFLO global-mesh analysis, are presented in Figures 18 through 24. Contours of the in- 
phase component or real part of the unsteady vorticity, source term and pressure are shown 
in Figures 18, 19 and 20, respectively, for vortical excitations at a — — x, —2 x and — 3x. The 
wave number magnitudes, |#c_oo|) and arguments with respect to the axial direction, a_oo, 
for these excitations are 5.08, 8.36 and 13.50 and —125.4 deg, —81.7 deg and —66.7 deg, 
respectively. As indicated in Figures 18 and 19, the vortical gusts are highly distorted as 
they are convected through the turbine blade row. The unsteady vortical and source term 
contour patterns for the gusts at a = — 2x and a = — 3x are quite different from those for the 
gust at a — —x. In the latter case the rectilinear vorticity contours far upstream of the blade 
row evolve into bowed shapes as the gust is carried through the blade row by the mean flow. 
The vorticity contours for the gusts at a = — 2x and a — — 3x are close to being straight lines 
far downstream of the blade row, but these show substantial changes in orientation relative 
to their upstream states. The source term contours in Figure 19 are severely distorted by 
the turbine blade row from mid-blade passage to the downstream boundary of the solution 
domain, particularly for the vortical gusts at a = — 2x and a — — 3x. Also, the source terms 
associated with the gusts at a = -2x and a = -3x have very large gradients within the 
blade passage and downstream of the blade row. These give rise to relatively large unsteady 
pressure responses, as indicated in Figure 20. Note that the vortical excitations at cr = — x 
and a = — 3x produce a subresonant pressure response; whereas the excitation at a = — 2x, 
a superresonant (1,1) response. 

The unsteady pressure- difference distributions acting on the reference blades of the tur- 
bine cascade, and its flat-plate counterpart, for vortical excitations at cr = — x, — 2x and 
— 3x are shown in Figures 21, 22 and 23, respectively. Moment responses to incident vortical 
gusts with vh,-oo • e*r = (1,0), u = 5 and -3x (-540 deg) < a < -x (-180 deg) are shown 
in Figure 24. The excitations at a = -396.8 deg and a = -311.7 deg produce resonant 
acoustic response disturbances far upstream of the turbine cascade, and far upstream and 
far downstream of the flat-plate cascade. Those at cr = —389.0 deg and at cr = —279.2 deg 
produce such response disturbances far downstream of the turbine cascade. The unsteady 
pressure differences associated with the turbine and flat-plate cascades are qualitatively sim- 
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ilar for a = — x, but differ substantially for a = — 2x and <7 = — 3x. The unsteady moment 
responses, particularly those for the in-phase moment, differ substantially over a significant 
range of interblade phase angles, i.e., —540 deg < a < —300 deg. Again, the flat-plate 
pressure-difference and moment responses predicted by LINFLO and the classical analysis 
of Smith are in very good agreement. 

The unsteady response predictions for the turbine and flat-plate cascades again indicate 
that mean-flow nonuniformities have a substantial impact on the unsteady aerodynamic re- 
sponse of cascades to vortical gusts, particularly for gusts occurring at high wave number. 
The large differences between the pressure- difference responses for the EGV and its cor- 
responding flat-plate cascade for a vortical excitation at a — — 3 x and between those for 
the turbine and its corresponding flat-plate configuration for excitations at a = — 2 x and 
a = — 3 x suggest that additional studies be undertaken to better understand the effects of 
mean-flow nonuniformity on blade response to high wave number vortical gusts. 

5.2 Viscous Flow: Acoustic Excitation 

Unstaggered Flat-Plate Cascade 

As a benchmark case to demonstrate the unsteady viscous-layer analysis, UNSVIS, we 
consider unsteady flow through an unstaggered flat-plate cascade. Here, the uniform inlet 
velocity is aligned with the mean positions of the blade chord lines. Hence, the steady 

background flow is uniform with velocity, pressure and total enthalpy given by V = e^, 
P = (7 A/ 2 ) -1 , and H — ~ + (7 — respectively. The unsteadiness is excited by 

an incident pressure disturbance, which carries acoustic energy toward the blade row from 
either far upstream (— 00 ) or far downstream (+ 00 ). The acoustic excitation travels in the 
axial (or chordwise) direction, i.e., k v = K y — a/G = 0 , and is therefore described by a 
velocity potential of the form [cf. (3.17)]. 

<£(x, t ) = -Re{[ 2 (/c etToo + w)] _ 1 p/.Too exp[i(/Cf Too £ -f ut))} , (5.3) 

where w, = k x and p/ iToo are the temporal frequency, axial wave number and complex 
amplitude of the unsteady pressure excitation, respectively. The axial wave number is related 
to the temporal frequency by [cf. (3.33)] 

= fu>M 

K t*°° ~ 1 ± m ’ 

Note that for a given temporal frequency, w, and Mach number, M < 1 , |/c^ +00 | > _<*,!, 

because a disturbance coming from downstream travels against the freestream flow direction. 

The flat-plate cascade produces no response to the pressure excitation described by 
Eq. (5.3); therefore, this equation describes the entire first- harmonic inviscid flow field. The 
time-dependent velocity, pressure, and total enthalpy in the flow are then given by 

P = P + p+ ... = ( 7 M 2 ) -1 + R e{p 1,+co expftX/c^^ + wt)] } + . . . , (5.5) 

V = V + v + ... = [l± MRe{p/ )=Foo exp[i(/c e>Too X + wf)]]}e { + . . . (5.6) 
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and 




H = i7 + /i + -- - =: -// + ( 1 i Jkf)R e {p/,=Foo Gxp[i(/c^^;ooC T <-><)]} + ♦ • ■ ■ (5*7) 

where p = -D^/Dt, v = V<£, and h = —D<f>/Dt + V<f> * are the first-harmonic unsteady 
pressure, velocity and total enthalpy, respectively. The real and imaginary components of 
the complex amplitude of the first-harmonic unsteady pressure acting along each blade and 
wake surface (0 ^ x = £ 2) for a cascade operating at M — 0.5 and subjected to upstream 

acoustic disturbances with pi,-oo “ (1 5 0) and cj = 1, 2, 5 and 10 are shown in Figure 25; 
similar results for downstream acoustic disturbances, in Figure 26. The pressure distributions 
shown in Figures 25 and 26 have been determined from Eqs. (5.4) and (5.5). 

The inviscid results in Eqs. (5.5) — (5.7) have been used as edge conditions for a series 
of unsteady viscous-layer calculations. The latter allow us to examine the effects of pressure- 
excitation amplitude, temporal frequency and axial wave number on the behavior of the 
viscous-layer displacement thickness, <5(x, i), surface shear stress f w = (Re) 1 (jJtdU / dn) w 
and wake centerline velocity Uqt . In particular, we can decompose each of these quantities 
into Fourier series; e.g., 

S(x,t) = / [1 -pU/(p e U e )]dn = Yj M x )exp(*W) , (5.8) 

^0 Ti~ — CO 

where 

/ j f 2 IT fw ^ 

<5„(x) = — — / <$(x, t) exp(—inu>t), n = 0,±l,±2, , (5.9) 

27r Jo 

to examine the behavior of their Fourier components. In this way we can gain insight into 
the relative importance of nonlinear viscous effects and, therefore, into whether a linearized 
viscous analysis could be applied to provide meaningful and efficient unsteady viscous-layer 
solutions. Note that the lower limit on the integral in (5.8) must be changed to -oo if the 
wake displacement thickness is to be determined. 

A series of calculations were performed for upstream and downstream acoustic excitations 
at the frequencies, a?, listed above and complex amplitudes, p/,± «>> of (0.1,0), (0.3,0), (0.5,0) 
and (0.75,0). In each case the steady Mach number is 0.5 and the Reynolds number, Re, 
is 10 6 . The viscous-layer calculations were initiated at (x,<) = (0.01,0), and the laminar 
flat-plate similarity solution was imposed as the initial condition in x and t. Also, the flows 
were assumed to undergo instantaneous transition from laminar to turbulent at x = 0.02. 

The viscous-layer calculation was carried out using 51 uniformly stretched (with = 
1.10) points across each boundary layer, with At/ = 0.0175 at the blade surface. The wake 
grid consisted of 101 points across the viscous layer stretched in the same manner as for the 
surface boundary layers. A total of 25 uniform time steps were used per temporal period of 
the unsteady excitation. Two different axial or streamwise mesh distributions were used — 
one for the excitations originating upstream and one for those originating downstream. In 
each case the streamwise distribution was selected so that there were at least 20 mesh points 
per wave length, 27 t//c € iToo , for the highest frequency considered, i.e., u = 10. 

For disturbances originating upstream, a variably spaced streamwise mesh was used with 
points clustered near the blade leading and trailing edges. The minimum streamwise spacing 
on the blade was Ax « 0.0177 at the blade edges, and the maximum was Ax « 0.0611 near 
mid-chord. The stretching used in the wake is identical to that used for the forward portion 
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of the blade, with Ax ~ 0.0177 in the first wake interval and monotonically increasing to 
Ax Ri 0.11 at the farthest downstream calculation point (x = 2.0). The resulting grid had 29 
points on the blade surface and 20 points along the wake. For the downstream disturbances 
a nearly uniform grid was employed, with Ax R* 0.025 on the blade and 0.025 < Ax < 0.030 
along the wake, where the grid was mildly stretched to distribute the points throughout the 
interval x € (1,2]. The resulting grid had 40 points along the blade and 36 points along the 
wake. 

First, we consider the solution for an incident pressure disturbance from upstream with 
Pi,- oo = (0.5,0), u = 5 and a = 0 deg. The temporal means, S 0 , r Wi0 and (7<g )0 , and the 
magnitudes, |<$ n |, l r w,n| and |f/<g )n |, of the first two harmonics (n = 1,2) of the displacement 
thickness, 6, surface shear stress, t w , and wake centerline velocity, Uq, as determined by 
the unsteady viscous-layer solution, are presented in Figure 27 along with the corresponding 
steady (|p/,-oo| = 0) results. The steady viscous solution has also been determined using the 
UNSVIS code, with the steady displacement thickness and surface shear stress given by 

~S(x) = /”[ 1 - pU/(hV.)]dn (5.10) 

J o 

and f w = (Ke)~ l (pdU /dn) w , where p and U are the density and streamwise velocity, respec- 
tively, in the steady background flow. The differences between the steady and the temporal 
mean values of the unsteady viscous quantities, and the amplitudes of the higher harmonic 
unsteady quantities provide a measure of the importance of nonlinear unsteady viscous ef- 
fects. 

The results in Figure 27 indicate that nonlinear effects are relatively unimportant in the 
unsteady flow driven by the prescribed upstream pressure excitation. However, similar results 
in Figure 28 for an acoustic disturbance originating downstream, with p/ +0 0 = (0.5, 0), ui = 5 
and <7 = 0, indicate that the nonlinear response of the viscous-layer displacement thickness 
can be quite significant. An unexpected result of the latter calculation is the predicted 
increase in the time-mean of the unsteady displacement thickness with increasing distance 
along the wake (i.e., as x — ► 2.0). In an attempt to determine whether this effect is physical 
or numerical in origin, an unsteady viscous solution was calculated using a grid with twice 
as many uniformly distributed streamwise points. It was found that, although slightly less 
pronounced, this behavior was also present in the fine grid calculation. 

The behavior of the zeroth- and first-harmonic components of the viscous parameters 
is illustrated in Figures 29 and 30 for upstream acoustic excitations with pi-oo — (0.5,0), 
(7 = 0, and uj — 1,2,5 and 10. Similar results for downstream acoustic excitations are given 
in Figures 31 and 32. For the upstream case, the time-mean values and the magnitudes 
of the first-harmonic viscous-layer response quantities are relatively unaffected by changes 
in excitation frequency. In addition, the first-harmonic displacement thickness is essentially 
out-of-phase with the first-order velocity fluctuation, u e , at the edge of the viscous layer, 
and the first-harmonic wall shear stress and wake centerline velocity are essentially in phase 
with the edge-velocity fluctuation. The corresponding results (Figures 31 and 32) for the 
downstream acoustic excitation show that frequency has a greater impact on the temporal 
mean and the first-harmonic amplitudes of the viscous quantities. In particular, the time- 
mean displacement thicknesses vary significantly with u and indicate that nonlinear effects 
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are important at higher frequencies. As indicated in Figure 32, the first-harmonic displace- 
ment thicknesses lag the fluctuating component of the viscous-layer edge velocity by phase 
angles ranging from approximately -140 to -180 deg; the first-harmonic wall shear stresses 
lead this fluctuating velocity by angles ranging from 0 to 40 deg; and the minimum wake- 
streamwise velocities are approximately in-phase with the wake-edge velocity fluctuation. 

As a final illustration, we consider the response of the viscous layer to different amplitudes 
of acoustic excitation. In particular, numerical results are presented in Figures 33 and 34 
for excitations at w = 5, a = 0 and p IiToo = (0.1,0), (0.3,0), (0.5,0) and (0.75,0). For 
disturbances from upstream the effect of unsteadiness on the time-mean values of S and 
t u is small, except for the highest disturbance amplitude, |p/,_oo| = 0.75, considered. The 
effect on the time-mean wake centerline velocity is negligible. The results for the acoustic 
excitations from downstream in Figure 34 once again indicate that such disturbances produce 
much stronger nonlinear viscous-layer responses, particularly in the displacement thickness. 
For both the upstream and downstream acoustic excitations, the phase angles (not shown) 
of the first-harmonic viscous quantities relative to the edge velocity fluctuation are almost 
independent of the excitation amplitude. 

The foregoing results indicate that the viscous-layer response parameters f u and Uq 
behave essentially in a linear manner for acoustic excitations originating upstream of the 
blade row. For excitations originating downstream, nonlinear effects can be significant at 
high temporal frequencies (say u> > 5) and/or high excitation amplitudes (|p/, +00 | > 0.5). 
For M = 0.5 the magnitude of the axial wave number for a pressure excitation from down- 
stream is three times that for an excitation from upstream. This produces a corresponding 
ratio in the magnitudes of the pressure gradients that drive the viscous solution, and is per- 
haps responsible for the nonlinear response of the viscous layer to downstream disturbances. 
Another factor is the direction of propagation of the acoustic excitation. We have performed 
a series of calculations to isolate these effects and have found the latter, i.e., direction of 
propagation, to be the dominant factor in promoting a nonlinear response in the viscous 
layer. 

The detected behavior of the displacement thickness in the far wake, i.e., increasing with 
increasing x, for high frequency or high amplitude downstream disturbances is a somewhat 
disturbing aspect of the present unsteady viscous solutions. A series of laminar calculations 
were performed, and the results showed a similar behavior. Thus, the increase in displace- 
ment thickness in the far wake (x — ► 2) cannot be attributed to the turbulence model used 
in the present viscous analysis. 


Turbine Cascade 

In order to demonstrate the coupled LINFLO/UNSVIS weak inviscid/viscid interaction 
solution capability, we consider an unsteady flow through the turbine cascade of § 5.1. This 
flow occurs at a Reynolds number, Re, of 5.0 x 10 5 , an inlet Mach number, JV/_ <*,, of 0.19 and 
an inlet flow angle, fl_oo of 45 deg. The unsteadiness is excited by a pressure disturbance 
from upstream, characterized by pi-oo = (0.35,0), = <r/G = 0 and u> = 1.0. The 

inviscid steady Mach number field and Mach number distribution along a blade surface are 
shown in Figure 16. For this example the unsteady pressure is essentially in-phase with 
the upstream acoustic excitation; the in-phase component of the unsteady pressure field is 
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depicted in Figure 35. The pressure excitation produces a superresonant (1,1) response. 
Therefore, the pressure contours in Figure 35 indicate a much higher pressure amplitude at 
the upstream boundary than that due to the acoustic excitation alone. The temporal mean 
and the upper and lower bounds (envelope) for the inviscid surface velocity are shown in 
Figure 36. The first-harmonic surface velocity, as determined by the input value of pi-oo, 
is of relatively small amplitude for this case. However, for pressure excitations at higher 
amplitude, the viscous-layer separates on the blade suction surface just upstream of the 
trailing edge, thereby precluding a continuation of the unsteady viscous calculation into the 

wake. 

Since the turbine blade has a blunt leading edge, the unsteady stagnation-point analysis, 
described in § 4.3, was applied at each time step to provide the upstream velocity profiles 
needed to advance the viscous-layer solutions along the blade suction and pressure surfaces. 
For this case, the unsteady stagnation point motion is confined to a small interval of length 
2|r/| 0.0037. A series of representative velocity profiles for the periodic flow within the 

stagnation region at four different times (i.e., t — tt /2,w,3tv /2 and 2 tt) are shown in Fig- 
ures 37 and 38. Here the instantaneous streamwise- velocity profiles are presented in the 
body-fixed frame of reference, and the abscissa on each plot refers to the location on the 
blade surface at which the velocity profile is determined. The interval shown in Figure 37 is 
centered about the mean location of the stagnation point; that in Figure 38, about the in- 
stantaneous location, i.e., the location at which U e = 0 in the body-fixed frame of reference. 
Note that the velocity profiles are shown over a much narrower interval, [-0.02|r/|, 0.02|r/|], 
in Figure 38, where the profiles indicate that reverse flow occurs in the immediate vicinity 
of the instantaneous stagnation point location. 

The viscous-layer calculation for the reference turbine blade and its wake was performed 
assuming that instantaneous transition from laminar to turbulent flow occurs at r/r te = 0.05 
on both the pressure and suction surfaces of the blade. Here, r is the distance along the blade 
surfaces measured from the mean leading-edge stagnation point location, and the subscript 
te refers to the values of r at the blade trailing edge. The grid used in this calculation had 
77 points along the blade surface and 54 points along the reference wake streamline. It is 
a stretched grid with At 0.002 at the farthest upstream point and At « 0.0001 at the 
trailing edge. The largest value of At on the blade, i.e., 0.052, occurs near midchord. The 
streamwise intervals grow aft of the trailing edge from At ~ 0.0001 to approximately 0.083 
one chord length downstream of the blade row. This distribution ensures that there are 
at least 20 points per streamwise wavelength of the unsteady excitation. The viscous-layer 
calculation was carried out using 71 uniformly stretched points across each boundary layer, 
with K v = 1.045 and A r) = 0.04 at the blade surface. The wake grid consisted of 141 points 
across the viscous layer stretched in the same manner as on the blade surface. A total of 40 
uniform time steps were used per temporal period of the unsteady excitation. 

Results of the unsteady viscous-layer calculation are shown in Figures 39 through 41. 
Temporal mean values and upper and lower bounds for the displacement thickness and wall 
shear stress along the upper and lower surfaces of the reference turbine blade are shown in 
Figure 39; corresponding results for the wake displacement thickness and minimum stream- 
wise velocity, in Figure 40. Here, the upper and lower bounds of a viscous-layer response 
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quantity, say the displacement thickness, are defined as 

~± oo 

6 = S 0 ± 2 ^2 l^„l , (5-11) 

71 = 1 

and 8w is the displacement thickness of the complete wake. Wake velocity profiles at four 
different instants of time are depicted in Figure 41. The unsteady response of the viscous 
layer is essentially linear for this example, i.e., the temporal mean and the steady viscous 
solutions are almost identical, and the Fourier amplitudes of the higher in > 2) harmonic 
components of <5, t w and t/w min are negligible. 

The foregoing results demonstrate the new capabilities that have been added to the 
UNSVIS code, i.e., an unsteady stagnation region analysis and a modified unsteady wake 
analysis. These results also demonstrate the present weak inviscid/ viscid interaction pre- 
diction capability that results from a sequential coupling of a linearized unsteady inviscid 
solution, determined using LINFLO, and a nonlinear unsteady viscous-layer solution. Un- 
fortunately, since boundary-layer separation usually occurs in realistic configurations, and a 
weak interaction analysis breaks down in such cases, a weak interaction analysis has a lim- 
ited range of application. This is particularly true for compressor cascades where, because 
of adverse mean pressure gradients, separation occurs near the blade leading edges. There 
is an important need, therefore, to develop a simultaneous coupling (or strong interaction) 
solution procedure for separated unsteady cascade flows. 
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6. Concluding Remarks 


Linearized inviscid and nonlinear viscous-layer unsteady aerodynamic analyses and codes 
have been developed for predicting subsonic unsteady flows through isolated, two-dimensional 
blade rows. The inviscid analysis can be applied to predict the aerodynamic response to 
structural (i.e., blade motions) and external aerodynamic (entropic, vortical, and acoustic) 
excitations. It provides very efficient response predictions for realistic cascade configura- 
tions, and should therefore be useful for turbomachinery aeroelastic and aeroacoustic design 
applications. The viscous analysis can be applied to predict the viscous-layer response that 
arises from imposed inviscid conditions at the blade and wake surfaces. At present only a 
weak or sequential coupling of the inviscid and viscous-layer solutions has been considered. 
The development of a strong inviscid/viscid interaction analysis involving a simultaneous 
coupling of the inviscid and viscous solutions is planned as future work. 

Under the present phase of this overall research effort a composite (global/local) mesh 
analysis has been developed for determining the linearized unsteady velocity potential. The 
composite- mesh capability allows a great deal of flexibility in prescribing a suitable mesh 
for the accurate resolution of an unsteady cascade flow, and this type of capability will 
be essential for analyzing gust encounters in unsteady transonic flows with moving shocks. 
However, based on the parametric studies reported herein, a global-mesh analysis appears 
to be quite adequate for analyzing low speed flows, provided that a sufficiently dense mesh 
is employed. Moreover, a global mesh analysis requires only about one- fourth of the CPU 
time needed for a composite analysis. Also, under the present effort, an existing unsteady 
viscous layer analysis has been extended by incorporating a similarity analysis for the flow 
in the vicinity of a moving stagnation point, and by properly accounting for the jumps in 
the inviscid flow variables across vortex-sheet unsteady wakes. 

The linearized inviscid (LINFLO) and nonlinear viscous-layer (UNSVIS) analyses have 
been demonstrated via application to compressor, turbine and flat-plate blade rows. The 
numerical results indicate that mean-flow nonuniformity has an important impact on the 
pressure responses of cascades subjected to vortical excitations. In particular, high wave 
number vortical gusts produce pressure responses in real blade configurations that differ 
substantially from those occurring in corresponding flat plate configurations. The viscous- 
layer analysis has been applied to unstaggered flat plate and to turbine cascades subjected 
to acoustic excitations coming from upstream or downstream and traveling in the axial di- 
rection. The flat-plate example is, perhaps, the simplest unsteady cascade problem that can 
be analyzed, both because of the geometric simplicity and because the unsteady pressure is 
nonsingular at the flat-plate leading edges. The numerical results indicate that the viscous 
layer responds linearly, for the most part, to acoustic excitations from upstream, but sig- 
nificant nonlinear response components occur for downstream excitations at high temporal 
frequency and/or high amplitude, which travel against the mainstream velocity. The numer- 
ical results for the turbine demonstrate the present weak inviscid/ viscid interaction solution 
capability for a realistic cascade configuration. 

Because of boundary-layer separation the range of application of a weak inviscid/viscid 
interaction analysis is severely limited. For example, the mean pressure rise produced by a 
compressor blade row typically causes boundary-layer separation very near the blade leading 
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edges, thereby precluding a continuation of a direct viscous-layer calculation along the blade 
and into the wake. Thus, the development of a strong inviscid/ viscid interaction analysis 
will be needed so that the effects of viscous displacement on the unsteady pressure response 
can be included in aeroelastic and aeroacoustic design calculations. The linearized inviscid 
response to high wave number gusts, and the linear vs nonlinear response of viscous layers 
to entropic, vortical and blade-motion excitations are other issues that will require further 
study and clarification. 
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List of Symbols 


All physical parameters listed below are dimensionless. Lengths have been scaled with re- 
spect to blade chord, time with respect to the ratio of blade chord to upstream freestream flow 
speed, density and viscosity with respect to their upstream freestream values, respectively, 
velocity with respect to the upstream freestream flow speed, normal and shear stresses with 
respect to the product of the upstream freestream density and the square of the upstream 
freestream flow speed, temperature with respect to the square of the upstream freestream 
speed divided by the fluid specific heat at constant pressure, and entropy with respect to 
the fluid specific heat at constant pressure. The number(s) in parentheses at the end of each 
symbol description indicates an equation in which the symbol appears. 

Roman 

«4woo Constant vector , (3.23). 

At Inner region turbulence model parameter, (4.11). 

A Speed of sound propagation in steady background flow, (3.7). 

a 0 Constant, (3.26). 

B Moving blade surface. 

B Blade mean position, (3.28). 

6, c Stagnation-region constants, (4.42), (4.49). 

dr Differential vector tangent to path of integration, (3.20). 

dry Differential element of arc length along a streamline, (3.19). 

e Unit vector, (3.18). 

T, V, Q Dependent variables in Levy-Lees transformation, (4.18). 

F Complex function of mean-flow stream function, (3.26). 

f,g Dependent similarity variables for stagnation region, (4.31). 

G Cascade gap vector (= Ge v ), Figure 1. 

H Total enthalpy, (3.48). 

H Total enthalpy in steady background flow, (3.37). 
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h 


First-order unsteady total enthalpy, (3.41). 
h! Turbulent total enthalpy fluctuation, (4.3). 

i Imaginary unit, (2.1). 

K{, Kr, Grid stretching factors used in viscous-layer calculations, (4.50), (4.51). 

/, / Functions appearing in viscous region similarity equations, (4.22), (4.34). 

M Mach number in steady background flow, (3.7). 

m Blade number index, Figure 1; complex amplitude of first-harmonic 

unsteady aerodynamic moment about midchord, Figure 14. 

n Unit normal vector, (3.3). 

P Pressure, (3.36). 

P Pressure in steady background flow, (3.7). 

Pr, Prj Prandtl number, turbulent Prandtl number, (4.3), (4.6). 

P+ Inner region turbulence model parameter, (4.11). 

p First-order unsteady pressure, Equation (3.40). 

p Complex amplitude of first-harmonic unsteady pressure, (2.4). 

Q Functions associated with Levy-Lees transformation, (4.17), (4.19). 

Q w Surface heat flux, (4.8). 

71 Surface (blade, wake or shock) displacement vector, (2.1). 

Re Reynolds number, (4.2). 

Re^ Reynolds number based on viscous-layer momentum thickness, (4.32). 

r Complex amplitude of surface (blade, wake or shock) displacement 

vector, (2.1). 

r T Flat-plate displacement tangential to surface, (4.30). 
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Sh Moving shock surface. 

Sh Shock mean position, (3.30) 

5 First-order unsteady entropy, (3.39). 

s Complex amplitude of first-harmonic unsteady entropy, (2.2). 

T Temperature, (4.4). 

T Temperature in steady background flow, (3.7). 

Tc Temperature constant used in Sutherland’s law, (4.5). 

t First-order unsteady temperature, (3.40). 

t Time, (2.1). 

U, V Unsteady velocity components along and normal to a moving blade 

surface or reference wake streamline, (4.1). 

Uq Wake centerline velocity, Figure 26. 

t/g Wake centerline velocity in steady background flow, Figure 26. 

f/<g >n Fourier component of unsteady wake centerline velocity, Figure 26. 

v' Streamwise and normal components of turbulent velocity fluctuation, (4.2). 

V Velocity, (3.1). 

V Velocity in steady background flow, (3.1). 

v First-order unsteady velocity, (3.4). 

v Complex amplitude of first-harmonic unsteady velocity, (3.5). 

Yr Complex amplitude of first-harmonic rotational unsteady velocity, (3.8). 

v. Complex amplitude of first-harmonic source-term unsteady velocity, (3.8). 

y' Turbulent velocity fluctuation, (4.14). 
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Cni< 


W 

W 

X 

x 

x,y 

y+ 

Greek 
ai,..., a 6 
0^00 
7 

7 ir 

A 

Ap 

A 6 A 77 


Moving wake surface. 

Wake mean position, (3.29). 

Lagrangian coordinate vector, (3.18). 

Position vector, (2.1). 

Cartesian coordinates along and normal to mean position of blade chord, 
Figure 1. 

Reynolds number based on friction velocity and normal distance from blade. 


Coefficient functions in Levy-Lees equations, (4.22). 

Acoustic attenuation constant, (2.4). 

Fluid specific heat ratio, (3.7). 

Turbulence intermittency factor, (4.15). 

Drift function, (3.19). 

Complex amplitude of first-harmonic unsteady pressure difference (lower 
surface — upper surface) acting on reference (m = 0) blade, Figure 7. 

Mesh spacings, (4.50), (4.51). 


6 

S n 

t 

e 

€ ii 

c 


Viscous-layer displacement thickness, (5.8). 

Viscous-layer displacement thickness in steady background flow, (5.11). 
Fourier component of viscous-layer displacement thickness, (5.8). 

Small parameter (<C 1); Turbulent eddy viscosity, (4.6). 

Turbulent eddy viscosity in steady background flow, (4.22). 

Turbulent eddy diffusivity, (4.6). 

First-order unsteady vorticity, (2.3). 
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c 

0 

0, A 

K 

P 

P 


f» V 


p 

p 

p 

a 

T 

T w 

T w 

T\u,n 
T, Tl 

$ 

P> 

<f> 


^«g,l ? ^ ^“e 


Complex amplitude of first-harmonic unsteady vorticity, (3.24). 

Cascade stagger angle, Figure 1. 

Dependent similarity variables, (4.32), (4.42). 

Wave number vector, (2.2). 

Molecular viscosity, (4.5). 

Molecular viscosity in steady background flow, (4.22). 

Cascade axial and “circumferential” Cartesian coordinates, Figure 1; 
independent similarity variables, (4.17), (4.33). 

Density, (4.1). 

Density in steady background flow, (3.6). 

First-order unsteady density, (3.40). 

Interblade phase angle, (3.5). 

Unit tangent vector, (3.3). 

Surface shear stress. 

Surface shear stress in steady background flow, Figure 26. 

Fourier component of surface shear stress, Figure 26. 

Viscous-layer coordinates along and normal to blade surface and 
reference wake streamline, (4.1). 

Velocity potential for inviscid steady background flow, (3.6). 

Velocity potential for first-order inviscid unsteady flow, (3.40). 

Complex amplitude of first-harmonic unsteady velocity potential, (3.11). 

Complex amplitude of first-harmonic convected potential, (3.25). 

Phase angles of the first-harmonic displacement thickness, wall shear stress, 
wake centerline velocity, and viscous-layer edge velocity, Figure 29. 
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X, Xo 


n 

u 

Mat hematical 

D/Dt 

D/Dt 

Im{ } 

Re{) 

V 

[] 

Subscripts 

B 

% 

E , R 

e 

I 

i, o 
/, n 
m 

N, T , z, 7 7 
and £ 


Outer region turbulence model parameters, (4.13), (4.14). 
Stream function for inviscid steady background flow, (3.20). 
Steady flow angle, Figure 1. 

Temporal frequency, (3.4). 


Convective derivative operator, (4.2). 

Convective derivative operator based on mean-flow velocity, (3.9), (3.40). 
Imaginary part of { }, Figure 7. 

Real part of { }, (2.1), Figure 7. 

Gradient operator, (3.6). 

Change in a flow quantity across a surface of discontinuity, (3.29). 


Reference blade surface, (2.1). 

Wake centerline, Figure 27. 

Acoustic excitation, response, (2.1). 

Edge of viscous layer, (4.7). 

Incident pressure or irrotational velocity disturbance, (2.4); 
initial r-station for viscous-layer calculation, (4.49). 

Inner, outer turbulent regions, (4.10), (4.12). 

Integer indices. 

Blade number index, Figure 1. 

Indicates direction: N — normal to inlet freestream direction, (3.18), 
T — parallel to inlet freestream direction, (3.18), z — out from the 
page, (3.20), t) — cascade “circumferential” direction, (3.31), and 
£ — cascade axial direction, (3.33). 
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5 

s 


te 


Moving blade, wake or shock surface, (3.2). 

Mean position of blade, wake or shock surface, (3.2); 
stagnation region variable, (4.28). 

Moving shock surface, (3.36). 

Blade trailing edge, (4.16). 


W, W min Wake value, (4.24); minimum wake value, Figure 39. 

w Wall, i.e., blade surface, value, (4.8). 

+ Axial locations at finite distance upstream and downstream from blade 

row, (2.4); point of intersection (£_,? 7 _) of axial line ( = and 
reference blade stagnation streamline, (3.20). 

ipoQ Far upstream/downstream freestream value of a steady flow variable, 

Figure 1; far upstream/downstream value of an unsteady flow variable, (2.4). 


Superscripts 

1 Turbulent fluctuation, (4.2). 

-(-j — Upper, lower viscous layer, (4.24). 

Steady background flow quantity, (3.11). 

Turbulent correlation, (4.2). 

— Time-dependent flow variable, (3.1), (3.4). 
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and <7 = 0: (a) t = 7r/2; (b) t = tt; (c) / = 3x/2; (d) t = 2%. 

Figure 38. Velocity profiles in the neighborhood of an instantaneous stagnation point 
location for the turbine cascade subjected to an upstream acoustic excitation with p/,-oo = 
(0.35, 0), uj = 1 and <7 = 0: (a) through (d) as in Figure 36. 

Figure 39. Temporal mean values and upper and lower bounds for the displacement thick- 
ness and surface shear stress along a turbine blade surface for an upstream pressure excitation 
with p/,_oo = (0.35,0), u = 1 and <7 = 0. 

Figure 40. Temporal mean values and upper and lower bounds for the displacement thick- 
ness and minimum wake streamwise velocity along a turbine wake for an upstream pressure 
excitation with p/,-oo = (0.35,0), u> = 1 and <7 = 0. 

Figure 41. Streamwise velocity profiles in the wake of a turbine blade for an acoustic 
excitation from upstream with pi-oo = (0.35,0), u> = 1 and <7 = 0. 
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Figure 2: Unsteady excitations: blade motion; incident vortical and entropic disturbances 
from upstream; and incident acoustic disturbances from upstream and downstream. 
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Figure 3: High Reynolds number steady flow over a blade surface. 
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Figure 4: Mach number contours and blade surface Mach number distributions for steady 
flow at M.oo = 0.3 and Q-oq — 40 deg through the EGV cascade. 
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Figure 5: Drift and stream contours for steady flow at M_oo 
through the EGV cascade. 


0.3 and Q_oo = 40 deg 
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Figure 6: Global and local meshes used for the EGV unsteady flow calculations. 
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Figure 7: In-phase component (real part) of the linearized unsteady potential due to the 
interaction of a vortical gust with v^-oo • e# = (1, 0), u = 5 and a = —2ir and the EGV 
cascade: (a) global-mesh calculation; (b) composite-mesh calculation. 
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Figure 8: Global- and composite- mesh solutions for the unsteady pressure-difference response 
due to the interaction of a vortical gust with VR t -oo • e# = (1,0), u> = 5 and a = — 2k and 
the EGV cascade: (a) in-phase component (real part) of the unsteady pressure difference; 
(b) out-of-phase component (imaginary part). 






Figure 9: Contours of the in-phase component of the unsteady vorticity for the EGV cascade 
subjected to vortical gusts with vr-co • e/v = ( 1 , 0) and u> = 5. 
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Figure 10: Contours of the in-phase component of the source term for the EGV cascade 
subjected to vortical gusts with 'Vr,- (X , • ©at = (1» 0) and u> = 5. 
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Figure 11: Contours of the in-phase component of the unsteady pressure for the EGV cascade 
subjected to vortical gusts with 'Vr^-oo ' = (1,0) and u) = 5. 
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Figure 12: Unsteady pressure-difference response of the EGV and corresponding fiat-plate 
cascades subjected to a vortical gust with vr-oo - ejv = (l,0),tu = 5 and a = -r. 
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Figure 13: Unsteady pressure-difference response of the EGV and corresponding flat-plat 
cascades subjected to a vortical gust with v H _ TO • ejv = (M)> u = 5 and ° = _2?r - 




Figure 14: Unsteady pressure- difference response of the EGV and corresponding flat-plate 
cascades subjected to a vortical gust with v R • e N = (1,0), u = 5 and a = - 37 t. 
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Figure 16: Mach number contours and blade surface Mach number distributions for steady 
flow at M_oo =0.19 and flLoo = 45 deg through the turbine cascade. 






Figure 17: Drift and stream contours for steady flow at M_oo — 0.19 and D.qo — 45 deg 
through the turbine cascade. 
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Figure 21: Unsteady pressure-difference response of the turbine and corresponding flat-plat 
cascades subjected to a vortical gust with 'Vr-oo • e/v = (1,0), w = 5 and a = — 7T. 




Figure 22: Unsteady pressure-difference response of the turbine and corresponding flat-plate 
cascades subjected to a vortical gust with _«> • — (1, 0), cj = 5 and <r = —2 it. 
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Figure 24: Unsteady moment vs interblade phase angle for the turbine and corresponding 
flat-plate cascades subjected to vortical gusts with Vr.-oq • ejv = (1 > 0) and u = 5. 
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Figure 26: Unsteady pressure distributions along blade and wake for an unstaggered flat- 
plate cascade (with fi = 0 = Odeg, M_ = 0.5 and G = 1) subjected to a unit-amplitude, 
|p/,+oo| = 1 , pressure excitation from downstream with k,,+<» = &G 1 = 0. 
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Figure 27: Temporal mean and Fourier magnitudes of the displacement thickness, <5 , wall 
shear stress, f w , and wake centerline velocity, Uq, for turbulent flow through an unstaggere 
flat-plate cascade (ft = 0 = Odeg, G = 1, M = 0.5 and Re = 10 6 ) subjected to an incident 
pressure disturbance from upstream with p/,_oo = (0-5,0), u = 5 and a - 0 deg. 







Figure 28: Temporal mean and Fourier magnitudes of the displacement thickness, 6, wall 
shear stress, t w , and wake centerline velocity, Ug, for turbulent flow through an unstaggered 
flat-plate cascade (Q = 0 = Odeg, G = 1, M = 0.5 and Re = 10 6 ) subjected to an incident 
pressure disturbance from downstream with p/ +00 = (0-5,0), u = 5 and a = Odeg. 
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Figure 29: Fourier amplitudes of 6, f w and Uq for an unstaggered flat-plate cascade (fi = 
0 = Odeg, G = 1, M = 0.5 and Re = 10 6 ) subjected to pressure excitations from upstream 

with pi _oo = (0.5, 0) and a = 0 deg. 
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Figure 32: Relative phase angles of first-harmonic displacement thickness, wall shear stress 
and wake centerline velocity for an unstaggered flat-plate cascade (ft = 0 = 0 deg, G = 1 
M = 0.5 and Re = 10 6 ) subjected to pressure excitations from downstream with p/,+oo = 
(0.5,0) and a = Odeg. 




Figure 33: Fourier amplitudes of <5, and Uq vs amplitude, |pj f _oo|> f° r an acoustic exci- 
tation from upstream, with u = 5 and <r = Odeg, interacting with an unstaggered flat-plate 
cascade (fi — 0 — Odeg, G ~ 1, M = 0.5 and Re = 10 6 ). 




Figure 34: Fourier amplitudes of 6, t w an 
tion from downstream with, u = 5 and a 
cascade (fl = 0 = 0 deg, G — 1, M = 0.5 


Uq vs amplitude, |p/,+oo|) for an acoustic excita- 
0 deg, interacting with an unstaggered flat-plate 


wid Re = 10 6 ). 
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Figure 35: Contours of the in-phase component of the unsteady pressure for the turbine 
cascade subjected to a pressure excitation from upstream with pi-oo — (0.35,0), u = 1 and 

(7 = 0 . 
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Figure 36: Temporal mean values and upper and lower bounds of the inviscid surface 
(viscous-layer edge) velocity for the turbine cascade subjected to an upstream pressure ex- 
citation with pi -oo = (0.35,0), w = 1 and <7 = 0. 
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Figure 38: Velocity profiles in the neighborhood of an instantaneous stagnation point location 
for the turbine cascade subjected to an upstream acoustic excitation with p /t _ = (0.35,0), 
= 1 and <r = 0: (a) through (d) as in Figure 36. 
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Figure 39: Temporal mean values and upper and lower bounds for the displacement thickness 
and surface shear stress along a turbine blade surface for an upstream pressure excitation 
with pi-oo = (0.35,0), u = 1 and a — 0. 






Figure 40: Temporal mean values and upper and lower bounds for the displacement thickness 
and minimum wake streamwise velocity along a turbine wake for an upstream pressure 
excitation with pi _oo = (0.35,0), u = 1 and a — 0. 







Figure 41: Streamwise velocity profiles in the wake of a turbine blade for an acoustic exci- 
tation from upstream with pi-oo = (0.35, 0), u — 1 and <7 = 0. 


97 



A. Details of the Viscous-Layer Solution Procedure 


The viscous-layer equations are solved numerically using an implicit finite- difference ap- 
proach, which is described in this Appendix. The finite-difference approximations used to 
discretize the governing equations and the quasilinearization applied to the resulting system 
of nonlinear equations is discussed in section A.l. The recursion relations and the associ- 
ated coefficients needed to solve the block- tridiagonal system of equations are given in A.2. 
Finally, the wake analysis is described in A. 3. 


A.l Finite-Difference Approximations 


The partial derivatives appearing in the governing viscous-layer equations are all first 
derivatives, with the exception of those accounting for viscous effects in the streamwise 
momentum and energy equations, which introduce second derivatives with respect to //. In 
the present analysis first-order accurate backward differences are used to approximate the £ 
and t derivatives. First and second partial derivatives with respect to // are approximated 
using central differences. The notation used below is defined in the illustrations of Figures 
A.l and A. 2. The subscript n is the mesh-point index in the //-direction, so that (•)„ refers 
to a quantity evaluated at the location, //„. 

First derivatives of the dependent variables are written as 


and 


where 


and 


(A.l) 
(A.2) 

(A.3) 

=6-6, 

A t — ^2 — ^1 > 

Af] n — 1 = Vn T]n—l 

Af) n = 77 n +i - Tj n . 


A.2 f n ^3 ,n 

MJ n * a* ’ 

<9A\ ^ A2 [T1 — Ai [tl 

dt I ~ At 

v / n 

fb4\ ^ A2,n+1 ~ A^n-l 
dv / „ ~ Ar} n -i + A// n 


Here, the subscript 1 refers to the current ^-station, previous time level, 2 refers to the current 
6station and current time level (i.e., the solution station) and 3 refers to the previous 6 
station, current time level (see Figure Al). Second derivatives in the //-direction are written 
in the form 



dr} 2 j n ~ A //„ -I- Aj/n-! \ 


^ (A n+ i A n ) (A n An—i)^ 
A 7? n A7/ n _! ) 


(A.4) 


If, as in all of the calculations performed here, | Aj/„/A?/ n _i — 1| <C 1, the //-derivatives remain 
formally second-order accurate. 
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Figure A.l: Computational molecule in the and t-directions. 
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Figure A.2: Computational grid in the 77-direction. 


Before discretizing the equations, the continuity, Eq. (4.19), streamwise momentum (4.20) 
and energy (4.21) equations, respectively, are rewritten in the following forms, 

A 1 ^ + ^ + A w F+A n =0 , 

d( dr] 

+ A-sT + A\T 2 + A h Q - ( l 


dt 


di ' dt] 


and 







(A.5) 

)" 

(A.6) 

m 

= 0 , (A. 7) 
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where the coefficients Ai, ... , An are given by [see (4.22)] 

A 1 =2(Q/Q, A 2 = 2(/(QU e ), A 3 = ai a 3 , 

A 4 = a 4 + (a 2 a 3 )/2 , As = — ( a 4 + a 3 ) , A 6 = oqae , 

A? = a 3 -\- (oc 2 ote)/2 , As = — (<*6 + <*s) , Ag = a 2 , 

A 10 = QjQ and An = — 0 * 10:3 . 

The nonlinear terms are quasilinearized by writing terms of the form AB as 


(A.8) 


AB « AB + AB - AB , 


(A.9) 


where the overbar denotes a known quantity (e.g., from the previous iteration). 

After quasilinearizing and discretizing the streamwise momentum equation, it can be 
written as 


where 


and 


^u,n^2.«— 1 + + ^n^2,n+l + B^ n V 2t „ + B° n Q 2 , n — Ai,» , 


bL 


- 2/ n-i 


'2,n 


Ai]n-i(Aj] n + A^.i) A T) n + Ai7„_i 

= ^7 + A(2^i,. - **.) + 4, + 2.4*?,, 


+ 


A^ n + A T), 


f ( n t l . 2 . ^ 

[Ar/ n Aij n _i J 


(A.10) 


(A.ll) 


(A.12) 



'2,n 

^‘n+1 

(A.13) 

Aljn 

Ar} n {Ai) n + Ai/„_i) ’ 

<„ 

_ F 2,n+l “ F \ 2,n—l 


(A. 14) 

Al)n + Ar/ n _x 

7 


= ^5 


(A.15) 

D„ 

^1^2, n -T, ^2,n+l - ^2,n-l , ^- 2 

= Ai + Ai + Vl '" A,„ + A,„_, + ^-- 

(A- 16) 


The energy equation is similar in form to the momentum equation, except that off- 
diagonal terms in both T and Q arise due to the presence of ^-derivatives of both quantities, 
resulting in the quasilinear discretized form 


^E,n^ • 2,n-l + B En T-. 2 ,„ + C En 7 r 2 ,n+ 1 + B En V 2>n + A E n C/ 2>n _i + B En G 2 , n + C En G 2 ,n+\ - D E<n 

(A.17) 


where 


.p _ 2Ag(l n _l / n -l)^2,n-l 

E ’ n - A7/ n _ 1 (Ar/ n + Ari n -i) ’ 

Be,ti = ^r(^2,n - £?3,n) + A 6 <?2,n + 2A 7 ^2,n 


(A.18) 
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2Ag(l n+ l ~ l n+ l)^2,n 2A 9 (/ n _l ~ l n ^x)^2,n 

Ar} n (Arj n + Arjn-i) ~ Arf^A^ + A r Jn . 1 ) ’ 


F _ 2A 9 (/ n+ l - / n .|.l)^ r 2,n+l 

B,n _ Arj n (Ar] n + Arj n _i) 

j-.y ^2,n+l £?2,n— 1 

E,n ~ Ar/ n + At ?„_! 


1 G - 
i E,n “ 


-V 2 ,n 


2/ n _i 


At;,, + At/„_i AT; n _ x (AT; n + At/„_i) 


flf, = A + A^, + ^, + 


2'n+t 


+ 


K-i 


c c , L = 


Aijn-iiArin + At/„-i) 
V 2 ,n 


Ar) n (Ar] n + Ar} n -i) 
+ -^8 > 


^n+j_ 


£,n A^„ + At/„_i At/„(A ?/„ + A^ n _ x ) 


and 


A 2 £l, n , _^1_ ^ 
At Af ' 


0*. = + 


V 2 ,n 


A77 n + A^ n _x 


((/2,n+l — £ 2 ,n-l) 


- , ^9 A 9 (/ n+ i ^1+5) , ^-2 -r-2 ^ 

+4^i,a, + *»*?.» + a ,„(a,„ + ^- n) 


a 

~ ^n-\) ( ^ t-2 ^ 

-*■ s (y 2,n — s 2,n— 1 / * 


Ar) n -i{Ar} n + Ar} n -i) 

The continuity equation, which is linear, has the discretized form 

^c,n^2,n-l + B E n T 2,„ + V 2 ,„ “ V 2>n -1 = D c<n 


where 


+ 


F _ MVi , ^ioAtv-i 

Ac ’ n ~ 2A( 

dF _ A? 

F^c.n ^c.n 


and 


D c<n = -"T^r - (^3,n + ^3,n-l) ~ A n AT} n -i 


(A.19) 

(A.20) 

(A.21) 

(A. 22) 

(A.23) 

(A.24) 


(A.25) 

(A. 26) 

(A.27) 
(A. 28) 

(A.29) 


A. 2 Recursion Relations 

To facilitate the solution of the above equations, the following recursion relations are 
introduced for Ti,n and Q% t n- 


^2 ,n 

= P E + g^V 2 ,„-l + R F n K 2,n-l + T F G2,n- 1 

(A. 30) 

Ql,n 

= Pn + Qn Vj,»-1 + + Tf^.n-l • 

(A.31) 
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Substituting these expressions into the discretized governing equations and collecting terms 
gives 


and 

where 


+ #£„&,» + JO’,, = d ; - 

+ W 3 E „V 2 , n = £ E - - Ag, , 


(A. 32) 


Kn 

pF i /iF pF 

K, 

pG . aOtF ^F 
~ x \i,n ' u u,n i n+l > 

Hi 

= Bl n + ClQl +1 , 

Hi 

— n _ /nrF pF 

— ^u,n u u,n r n+l 


(A.33) 


and 


(A.34) 


i jE nF , pF pF , />G pF 

^ l,n — n E t n + L/£; n it n+1 + n /t n+1 , 

tjE pG , rfF t-iF . nG cpE 1 

n 2,n n E,n T U F,n J n+l T L 'E t n 1 n + 1 > 

Ht = ^ in + Cf tn Q^ 1+ C7f, n gf +1 , 

= ^,n - Cl n P* +x - Cf in P n E +1 . 

Solving for V 2]T ,, using continuity (A. 26), and substituting into the momentum and energy 
equations yields 

hC ? 2 .„ + # 2 *,X« = n: F - A;V 2 , n _! - Afv 2 , n _x - 

and 


respectively, where 


and 


= c;' 

- A*; F 2 ,n-l - A'J v 2 , n 

< 

ttF tjF pF 

~~ -“1 ,n ” ^3,n^c,n 5 

< 

= Hi, 

d-: 

= Dl-H[ n D c>n , 

< 

= ^n, 

F 

K 

— 4 ^ _ 4 ^ 

K 

= 0 

< 

r/F ttE pF 

“ ^M,n ^3,n^c,n > 

K 

pF 

— ^*2,n > 

d-: 

= d*-h°„d„, 

A’e 

pF 

” ^3 t n > 

a ’e 

— A F — h e a f 

™E,n ,n^c,n 5 

A'/ 

= ^E,n ' 


(A. 35) 
(A. 36) 


(A.37) 


(A. 38) 
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Letting bp and be represent the right-hand sides of the momentum and energy equations, 
respectively, leads to the system of equations 


1 ,» 

h£ ' 

1 .n 

K J 


(£ 1 -U)' 


which can be solved using Cramer’s rule. Defining A by 

a tj* f if* E tj* e h* f 

i.e., the determinant of the coefficient matrix, T 2 , n can be written in the form 


(A. 39) 


(A.40) 


= % { d ? - Afv,,-. - 


H: 


* (D-‘ - A'e V,,_, - A '/ - Ajfe,.-,) • (A.41) 


Comparing Eq. (A.41) to Eq. (A. 30) gives 



Pi = “nDf + a u D‘‘ , 
Qn — ~ a \\A* u — a 12 A E , 
K = -<hiA' u F - a 12 A* E F 

(A.42) 

and 

Tn = — a u A* u — a\ 2 A* E , 


where 

h£ 

a„ = A 

(A.43) 

and 

h£ 

012 = " A • 

(A.44) 


Analogous expressions are easily derived for Q 2 ,n- 

The recursion coefficents at the ^-location 7? n , i.e., Pf, P„, etc., are related to the values of 
the coefficients at points m > n only, while the expressions for the flow variables themselves 
relate the unknown variables at n to the known values at points m < n. Therefore, the 
recursion coefficients can be evaluated by starting at the outer edge, n — ne, using the edge 
boundary conditions to set the values of the coefficients there, i.e., 


PL 

= i, 

pE 

r ne 

= 1, 

Qle 

= o, 

Qt 

= o, 

RL 

= 0, 

Rt 

= 0, 

rpF 
1 ne 

= 0, 

rpE 
1 ne 

= 0, 


(A.45) 


and sweeping inward to the wall, solving for all of the coefficients. The flow variables are 
next determined by using the wall boundary conditions [i.e., Eq. (4.26)] to set T, V and Q 
at n = 1, and then directly applying the recursion relations to solve for the flow variables by 
sweeping outward from the wall to the edge of the viscous layer. 
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Figure A. 3 : Computational grid in the 77-direction for the wake. 

A. 3 Wake Solution Procedure 

The wake viscous-layer analysis is similar to the blade-surface analysis, with changes to 
the boundary conditions and the introduction of modifications to the definitions of some of 
the independent and dependent variables, as described below. Edge conditions are applied 
at both the upper and lower edges of the wake viscous layer and the 77 = 0 coordinate line 
is assumed to correspond to the wake reference streamline, so that V = 0 at 77 = 0 . As 
discussed in Chapter 4 , the variables used in the upper- (i.e., 77 > 0 ) and lower-wake (i.e., 
77 < 0) regions are defined differently, so that, for example, in the upper wake 

T = T* = and Q = Q + = -£- , (A. 46 ) 

U+ H+ 

while in the lower wake 

f = f- = S- and Q = Q- = -§~. (A. 47 ) 

U e H c 

In order to ensure that the physical variables (i.e., U, V and H) remain continuous across 
the entire wake, the different variable definitions must be accounted for at the interface 

between the upper and lower portions of the wake, at 77 = 0 . This is easily done when 

evaluating the governing equations at 77 = 0, where both upper- and lower-wake variables 
are used simultaneously. 

The details of the wake solution procedure are given below. For the purpose of the 
present discussion, we assume that the mesh-point index in the 77 direction, n, is equal to 
one at 77 = 0 and increases in the 4-77-direction (i.e., in the upper wake), with n = ne at the 
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upper edge of the 77-mesh. In the lower wake, the mesh-point index is given by 21, 31 ,. . ., 
nel, ascending with increasing distance from the reference wake streamline (see Figure A.3), 
with n = nel at the lower edge of the 77-mesh.. 

Wake Boundary Conditions 

The boundary conditions for the wake are given in Eqs. (4.25) and (4.27). The imple- 
mentation of these conditions into the numerical algorithm is discussed below. 

The recursion relations used in the wake are defined like on the blade surface, i.e., in the 
upper-wake region, 


and 


^2,» — Pn + Qn^2,n-1 + ^2,71-1 + T,fV/2,n-l 

Qi ,n = Pn + QnV"i,n-l + P-n'P- 2,n-l + > 


and in the lower- wake region, 

■T 7 2 ,nl = Pnl + Qni^2,{n-\)l + P-nl^ + ^rU^2.(n-l)/ 

an ^ @2,nl ~ Pnl + QnlVl,{n-l)l + Rnl^‘. 2 , (»-!)/ + Tnl^2,(n-l)l • 


(A. 48) 
(A.49) 


(A. 50) 
(A.51) 


With the edge conditions on T and Q given at n = ne and n = nel, the edge values of the 
recursion relations are written as 


PL = PL = P* = PL = 1 , ( A .52) 

Qi = Qi = C = Ri = TL = Tg = 0 (A.53) 

and 

Qi, = Qi, = Ri, = Ri, = TL = tL -- 0 . (A.54) 

The recursion relation coefficients are determined by applying Eqs. (A.42) in both the upper- 
and lower- wake regions. 

The values of T and Q at the wake reference streamline (n = 1) are determined by 
writing the momentum and energy equations there, substituting the recursion relations from 
the points immediately above and below this line (i.e., at n = 2 and n — 21, respectively), 
and simultaneously solving the resulting equations. The governing equations are written at 
n = 1 in terms of upper-wake variables, so that whenever a value of T or Q at n — 1 = 21 is 
used, it is transformed to the appropriate form. Thus, at n = 1, 


P 7 -. 2,n— 1 = JTtP 7 . 2 , 2 1 an d £/2,n-l = -~~rQ2,2l • 

VJ tij 


(A. 55) 


The discretized s-momentum equation at the point n = 1, with the boundary condition 
V24 = 0 incorporated, is 


+ < 1^.1 + + K-fiv = D -.i • 


(A.56) 
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The recursion relations for T , written at n = 2 and n — 21 after expressing T and Q at n = 1 
in terms of upper variables, are given by 


and 


2,2 — P 2 + P-l? 2,1 + 


^2,21 — P[l + R 2 ljfz ^2, 1 + 2,1 


(A.57) 


Using Eq. (A. 58) for ^ 2,2 in the momentum equation gives 



+ H[, 52,, = D? , - , 

(A. 58) 

where 

+ CjX , 

fl£. = B?, + C[,,Tf 

(A. 59) 

and 

Ol = A... - C F /' F . 


Similarly, substituting for ^ 2 , 2 ; in Eq. (A. 58) yields 



i/rV 2 ,i + ^ 2 *^ 2, 1 = , 

(A. 60) 

where 

#r F =^ + < 1^7 




(A.6i) 

and 

£>;' - j - A f , ? e . 

’ ’ t/u+ 


The 

same procedure is applied to the energy equation, which is written at n 

= 1 in the 


form 


2,2/ + Bej'P- 2,1 + C El T-. 2,2 + Agj -*+@2,21 + ^£,1^2,1 + C El G2,2 — -Df;,! 5 (A.62) 


where V 2 ,i has been eliminated. The recursion relations for Q, written at n = 2 and n = 21 
after expressing T and Q at n = 1 in terms of upper variables, are given by 


and 


£ 2,2 


@ 2,21 


= J P F + ^2, 1 +T ; f^,l 

77 + 

jjE p£ 77 1 nnE 

= “ 2 J + 2,1 + ^2/ 



^2,1 • 


(A. 63) 
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Substituting for J- 2,2 and £ 2,2 in Eq. (A. 62) yields 


where 


and 


Hi, 1^2,1 + Hf \Q 2 ,1 — Df — Ag j ' A ^2,11 ~ A% j ^+$2,21 , 


Hi 1 — Bf t + Cf t Rf + C El Rf , 

ttE tdG 1 riF rpF . s~iG rpE 

M 2,1 “ ^£,1 + L / E , l 1 2 ~ T K -' E ,\ 1 2 

Df = D ea - C^Pf - C% x Pf . 


(A.64) 


(A. 65) 


Substituting for ^ 2,21 and Q 2 , 2 i gives 




Hf T 2,1 + H?Q %X = Df 

(A. 66) 

where 

Hf 

- Hi,i + A E iR 2i + A E1 R 2l ~^_ , 



H? 

riE , A F rpF^e H e .q j, E 

- n 2l + a E 1 i 2l ^ + a e i i 21 

( A.67) 

and 

Df 

- n E A F P F EL A a pE He 

~ U \ ~ A Ey2l jj+ A E,\ r 2l jj+ ■ 



With the starred coefficients defined in this form, Eqs. (A.60) and (A. 66) can be solved as 
before, to yield 


and 


where 


T%\ 

£ 2,1 


H±_ _ n * F _ 

A Dx A ° x 


Ell D .r , E? 

A 1 + A 1 ’ 


a = h; h* -h; h* 


(A.68) 

(A.69) 

(A.70) 


Note that these reference streamline values are expressed in terms of upper-wake variables. 

The solution in the wake is computed as follows. Using the edge conditions for the upper- 
and lower-wake viscous layers, the recursion relation coefficients for each are computed by 
sweeping inward from the edges to the centerline. The flow variables T and C? 2 ,i are 
calculated from Eqs. (A.68) and (A.69), yielding their upper-wake values (i.e., and 
$ 21 )- The quantities T and Q at n = 1 need to be re-expressed in in terms of the lower- wake 
variables, in order to apply the recursion relations, Eqs. (A. 50) and (A.51), to obtain the 
solution in the lower-wake region. These are calculated from the expressions 




and Q 2l = JflGt 1 


(A.71) 
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With the variables at n = 1 known, including Vi t \ from the boundary condition applied 
there, the recursion relations and the continuity equation are used to solve for all of the 
unknowns at each mesh point, sweeping from n = 2 to n = ne on the upper wake and from 
n = 2/ to n = nel on the lower wake. 

Because the nonlinear discretized equations have been quasilinearized, the solution pro- 
cedure described above is applied iteratively at each spatial and temporal location, until the 
local solution has converged to a specified tolerance level. 
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